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ABSTRACT 


Yucuis,  William  Allen.  M.S.,  Purdue  University,  December  1982. 

Computer  Simulation  of  a  Multiaxis  Air-to-Air  Tracking  Task  Using 
the  Optimal  Pilot  Control  Model.  Major  Professor:  David  K.  Schmidt. 

The  primary  objective  of  the  research  is  to  simulate  the  multiaxis 
air-to-air  tracking  task  and  determine  if  the  augmentation  control 
synthesis  results  in  improved  piloted  vehicle  performance.  Also,  the 
simulation  results  will  be  compared  to  previously  obtained  analytical 
results  to  see  if  the  optimal  pilot  control  model  should  be  modified. 
The  simulation  was  accomplished  on  the  Aeronautical  Engineering  mini¬ 
computer  simulation  system  and  the  results  substatiate  the  claim  of 
improved  performance.  Since  the  analytical  model  and  simulation  task 
have  different  target  motions,  nothing  definitive  can  be  said  about 
modifying  the  model. 
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CHAPTER  1 .  INTRODUCTION 


1.1  Motivation 

The  trend  in  future  fighter- type  aircraft  design  is  towards  the 
designing  and  building  of  lighter,  more  maneuverable  aircraft.  To  attain 
these  goals,  designers  have  turned  to  new  advanced  in  technology,  such 
as  active  flutter  suppression  control  systems,  direct-lift  or  side-lift 
control  surfaces,  and  lighter,  yet  stronger  structural  materials.  The 
use  of  such  advanced  methods,  however,  radically  alters  the  performance 
characteristics  of  the  piloted  vehicle,  and  classical  methods  for  pre¬ 
dicting  these  characteristics  have  proved  inadequate.  The  classical 
methods  use  some  type  of  specification  on  the  vehicle  rigid  body  mode 
that  has  traditionally  resulted  in  acceptable  piloted  vehicle  performance. 
The  advanced  methods  mentioned  above,  however,  introduce  higher  order 
dynamics  into  the  system,  which  are  assumed  insignificant  in  classical 
approaches,  and  are  a  cause  for  errors.  Therefore,  some  other  method, 
or  combi natij.i  of  methods,  must  be  used  to  accurately  predict  the  piloted 
vehicle  performance  when  higher  order  dynamics  are  present.  One  such 
approach,  which  will  be  tested  in  this  research,  uses  the  optimal  pilot 
control  model  (OCM),^  to  directly  predict  piloted  vehicle  performance. 
The  main  purpose  of  this  research  is  to  simulate  a  multiaxis  air-to-air 
tracking  task,  which  includes  higher  order  dynamics,  and  to  compare  the 
simulation  results  to  already  available  analytical  predictions.  The 
aim  Is  to  validate  the  control  synthesis  methodology  used  to  obtain  an 


augmentated  control  system  and  to  determine  if  any  changes  to  the  OCM 
are  required.  The  analytical  development  of  the  control  synthesis  has 
been  done  by  Schmidt,^ and  provides  the  theoretical  foundation  for 
the  research. 

1.2  Outline 

The  report  will  begin  by  presenting  an  analytical  development  of 
the  optimal  control  pilot  model  and  then  explaining  the  control  law 
synthesis  used  to  develop  an  augmentation  system.  Chapter  2  will  review 
the  basic  OCM  and  then  apply  the  control  synthesis  to  a  pitch  tracking 
and  a  multiaxis  tracking  task.  The  analytical  study  involves  a  fighter 
type  aircraft  with  two  different  sight  configurations,  which  introduce 
higher  order  dynamics  into  the  system.  Finally,  some  additional  target 
motion  is  added  to  the  basic  model  to  more  realistically  simulate  the 
task. 

Chapter  3  outlines  the  simulation  system  used  to  accomplish  the 
simulation.  The  various  hardware  components  and  their  capabilities  will 
first  be  present  and  then  the  development  of  the  computer  software  will 
be  discussed. 

Chapter  4  explains  how  the  multiaxis  tracking  task  simulation  is 
used  to  test  the  research  hypotheses.  The  assumptions  used  in  the 
simulation  are  examined  and  the  evaluation  process  is  explained.  An 
outline  for  operating  the  simulation  system  is  given,  along  with  the 
methodology  for  analyzing  the  data.  Finally,  a  summary  of  simulation 
and  analytical  results  is  tabulated  and  discussed. 

The  final  chapter.  Chapter  5,  uses  the  results  given  in  Chapter  4 
to  derive  some  conclusions.  The  results  validate  the  primary  hypothesis 
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of  improved  piloted  vehicle  performance  using  the  augmentation  control 
synthesis.  The  results  do  not  yield  any  definitive  conclusions  about 
changing  any  of  the  OCM  parameters,  since  the  addition  of  the  target 
motion  makes  an  accurate  comparison  impossible.  However,  the  simula¬ 
tion  results  do  give  some  promising  trends,  indicating  the  OCM  is  a 
good  model  of  the  task.  Finally,  some  suggestions  for  future  research 
are  given. 
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CHAPTER  2.  MODEL  DEVELOPMENT 

2.1  Introduction 

(51 

The  optimal  control  pilot  model  (OCM)'  '  has  been  used  successfully 

to  accurately  predict  piloted  vehicle  performance.^  Since  it  can  be 

readily  adapted  to  a  problem  that  includes  higher  order  dynamics,  the 

OCM  was  chosen  to  model  the  multiaxis  air-to-air  tracking  task.  An 

augmentation  control  synthesis  methodology  has  been  developed  by 
(21 

Schmidt  which  uses  the  optimal  pilot  control  model  as  the  system 
model.  This  chapter  will  present  a  brief  overview  of  the  OCM  and  the 
control  synthesis  and  then  apply  them  to  both  a  pitch  tracking  and  a 
multiaxis  tracking  task.  A  comparison  between  analytical  and  simulation 
results  for  the  unaugmented  pitch  tracking  task,  as  well  as  augmented 
versus  unaugmented  analytical  results  for  both  tasks,  will  show  the 
advantage  of  using  the  OCM  and  also  indicate  that  the  augmented  system 
could  have  significant  improvements  in  piloted  vehicle  performance. 
Finally,  some  changes  to  the  OCM  will  be  made  to  more  realistically 
simulate  the  target  vehicle  motion. 

2.2  Optimal  Pilot  Control  Model  and  Control  Synthesis 

The  multivariable  optimal  control  pilot  model  was  selected  because 
it  has  several  important  advantages.  First,  the  model  easily  adapts 
Itself  to  the  higher  order  dynamics  present  In  the  problem.  Secondly, 
the  OCM  optimizes  the  piloted  vehicle  performance,  which  is  the  ultimate 


goal  of  ar\y  aircraft  designer.  The  OCM  accomplishes  this  by  not  only 
including  the  vehicle  characteristics  in  the  model,  but  also  incorporat¬ 
ing  the  pilot's  "strategy"  for  achieving  the  desired  performance  criteria. 
Finally,  the  interactive  structure  between  the  pilot  and  the  augmenta¬ 
tion  control  system  is  embedded  within  the  model.  This  structure 
accruately  reflects  both  the  cooperative  nature  of  a  pilot  controlling 
an  augmented  vehicle  and  the  pilot's  ability  to  readily  adapt  to  various 
control  levels. 

The  optimal  control  pilot  model  represents  the  system  in  the  famil¬ 
iar  linear  form  of  x  =  Ax  +  Bup.  The  OCM  then  assumes  that  a  highly 
trained,  well  motivated  pilot  will  adopt  a  "strategy"  of  choosing  the 
control  inputs,  up,  to  minimize  an  objective  cost  function,  Jp,  which 
is  given  by 

Jp  «  E{lim  }  jV Qy  +  upRup  +  ujGJp)dt}  (2.2.1) 

The  matrices,  Q  and  R,  are  weightings  on  the  pilot's  observation  vector, 
y,  and  the  control  inputs,  up,  while  the  matrix  G  is  selected  to  obtain 
a  specified  neuromuscular  lag,  tN.  ^  A  control  law  synthesis  method 
is  used  to  improve  the  piloted  vehicle's  performance. The  aug¬ 
mented  system  is  first  represented  by  x  =  Ax  +  Bup  +  Bua  +  w,  where  ufl 
is  the  equivalent  control  system  input,  and  w  is  system  noise.  The 
augmented  control  inputs,  u_,  are  then  selected  to  minimize  an  augmented 

fl 

objective  cost  function,  J_,  which  is  given  by 

fl 

VV  Et’1m  T  (5I  F5.)dtl 


(2.2.2) 
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The  values  of  u.  that  minimize  J.  are 

fl  Q 

ua  =  -Kx  x  -  Ky  iip  (2.2.3) 

where 

Kx  =  F_1BTKA1  (2.2.4) 

Ku  '  F"1bT|(A2  (2.2.5) 


and  KA1  and  are  obtained  by  simultaneously  solving  two  coupled 
(21 

Ricatti  equations.  The  resulting  system  is  therefore  represented  by 


where 


x  *  W 


+  B  u 
aug  p 


A 

aug 

B 

aug 


A_BKx 

B*BKu 


(2.2.6) 


(2.2.7) 

(2.2.8) 


The  next  step  is  to  determine  the  unique  objective  cost  function  for 
the  particular  task  in  question  and  then  apply  the  augmentation  synthe¬ 
sis.  The  first  task  to  be  considered  is  a  pitch-only  tracking  task. 


2.3  Pitch  Tracking  Task 

The  initial  task  to  be  examined  is  a  pitch-only  tracking  task,^ 
where  a  fighter  type  aircraft,  employing  two  significantly  different 
tracking  sights,  attempts  to  track  a  wings  level  target.  This  task  was 
selected  because  the  Inclusion  of  the  sight  dynamics  introduce  higher 
order  dynamics  into  the  system.  Therefore,  with  the  same  vehicle  rigid 
body  dynamics,  the  introduction  of  higher  order  system  dynamics  will 
significantly  alter  the  piloted  vehicle  performance.  This  section  will 
first  explain  the  differences  between  the  two  sight  configurations,  and 
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how  the  higher  order  dynamics  are  Introduced  Into  the  problem,  and  then 
present  the  basic  and  augmented  model  parameters,  with  an  explanation 
of  how  the  parameters  were  selected. 

Any  aircraft  tracking  sight  system  uses  the  known  ballistics  of 
the  attacking  aircraft's  weapons  system  plus  the  available  data  for  the 
attacker  and  target  states,  to  "predict"  where  the  next  bullet  would  hit 
if  the  gun  were  fired  at  this  moment.  Therefore,  if  the  target  is 
located  within  the  sight,  and  the  gun  is  fired,  a  bullet  would  theoret¬ 
ically  impact  the  target.  The  accuracy  of  a  sight  system  is  therefore 
dependent  on  how  accurately  the  ballistics  and  states  are  known.  A 
typical  assumption  is  that  the  ballistics  and  the  attacker's  states  are 
correctly  known,  so  that  the  target  states  are  the  only  variables.  The 
two  sight  configurations,  then  represent  different  methods  of  calculat¬ 
ing  the  target's  states.  The  first  sight,  called  an  ideal  sight, 
assumes  the  target's  true  line  of  sight  rate, a,  and  its  true  normal 
acceleration,  aT,  are  known.  The  displayed  sight  lead  angle,  x,  is 
then  calculated  by 

T  ^  yy 

X  =  Tf(©-6)  -  2§-  aT  -  -038  a  (2.3.1) 

where  e,  a.  and  V  are  the  attacker's  pitch  rate,  angle  of  attack,  and 
velocity,  respectively,  D  is  the  target  range,  and  Ty  is  the  projectile 
time  of  flight.  ^  The  second  sight,  called  the  typical  sight,  approxi¬ 
mates  e  and  aj  by  assuming  g  equals  x»  the  displayed  lead  angle  rate, 
and  aT  *  Z^o,  the  attacker's  normal  acceleration.  The  pitch  tracking 

task  involves  minimizing  e,  the  tracking  error,  where  c  *  x-g.  using 
2 

a  simple  K/s  plant.  The  pitch  tracking  task  display  is  shown  in 


Figure  1.  The  fixed  reference  point  represents  the  weapon's  boreline, 
while  the  sight  is  dynamic,  and  moves  in  the  display.  It  is  these  sight 
dynamics  that  introduce  the  higher  order  systems  dynamics  into  the 
problem.  By  comparing  the  transfer  functions  of  e(s)/6st(s)  for  the 
two  sight  configurations,  as  given  in  Table  1,  it  should  be  noted  that 
the  numerator  coefficients  are  functions  of  both  the  vehicle  and  sight 
characteristics,  while  the  denominator  coefficients,  except  for  Tf  in 
the  typical  sight  configuration,  are  functions  of  only  the  vehicle 
characteristics.  Therefore,  with  the  two  sights  introducing  higher 
order  system  dynamics,  a  unique  objective  cost  function  must  be  found 
that  allows  the  OCM  to  represent  this  system. 

The  problem  of  determining  a  unique  objective  cost  function  for  a 

(3\ 

pitch  tracking  task  was  done  by  Harvey.  '  The  pitch  tracking  task  was 
flown  using  several  actuator  dynamics  and  the  two  sight  configurations 
mentioned  above.  Harvey  then  used  fixed  base  simulation  data  to  deter¬ 
mine  a  unique  objective  cost  function.  Harvey's  simulation  data  was 
for  an  aircraft  with  the  typical  sight  configuration,  and  a  comparison 
between  analytic  and  simulation  results  is  presented  in  Table  2.  Note 
the  excellent  correlation  between  the  experimental  results  and  the  typi¬ 
cal  sight  analytical  predictions,  plus  the  improvements  resulting  from 
the  ideal  sight  configuration.  He  achieved  excellent  correlation  over 
a  wide  range  of  dynamics,  so  his  OCM  parameters,  as  shown  in  Table  3, 
were  used  for  applying  the  augmentation  strategy.  The  augmentation 
synthesis  will  now  be  applied  to  the  OCM  to  analytically  predict  improve 
ment  in  piloted  vehicle  performance. 

The  augmentation  control  synthesis  outlined  In  Section  2.2  was 
applied  to  both  the  ideal  and  typical  sight  configurations.  The 
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Table  1 

Transfer  Function  Comparison 

Second  order  plant: 


IT 


11.7 


Ideal  sight  display: 

sis)  ,  K<v3  *  *»;  v  *  V 

sst's|  S2(V  ♦  1)(sZ  ♦  ♦  » 


sp 


? 


Typical  sight  display: 

& 


KfMjS3  +  MgS2  +  MjS  +  M4) 


sz(Tfs  +  l)(xas  ♦  l)(sZ  +  2cspwsps  +  »sp>) 


K  *  -16  r/D 

a 

^  «  -.9627  Z6Tf 

N2  -  (D/V  ♦  Tf)(VMj  ♦  26M-)  -  Zfi 

N3  -  VM q  +  Ma>  +  ^  +  V(MaZ*  *  W 

N4  “  "aZ6  -  ZaM6 

M1  -  Z6Tf  (ZaO/2V2  -  .962) 

^  -  N2  +  lf{laD/Z^  -  .962)  (VMfi-  MqZ6) 


i 

i 

i 

i 


Table  2 


Table  3 


Pilot  Model  Parameters 

Observation  delay  ,  t  *  0.2  sec. 

Neuromuscular  time  constant,  =  0.1  sec. 

Observation  vector,  y  T  =  (e,  e,  x,  x) 

Cost  function  weightings,  Q  =  (16,  1,  0,  4) 

yii 

Ru  =  0 

Full  attentional  allocation 
Observation  thresholds,  T  =  T  =  0.65  deg 

£  A 

T*  =  T*  =  1.3  deg/sec 

£  A 

Observation  noise  variance,  V  =  irp„T^„ 

yi  y  yi 
Py  =  0.01 
2 

Motor  noise  variance,  Vu  =  hpuTu 

P  P 

pu  =  0.001 


calculated  gains,  plus  performance  results,  are  tabulated  in  Tables 
4  and  5,  and  as  hoped,  the  augmented  system  resulted  in  a  lower  error 
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e,  as  well  as  a  reduced  stick  deflection,  6st,  and  stick  rate,  6st> 
Therefore,  the  control  synthesis  predicts  not  only  an  improvement  in 
performance,  as  evidenced  by  reduced  error,  but  also  indicates  a  lower¬ 
ing  of  pilot  workload,  as  shown  by  reduced  stick  deflections  and  rates. 
The  control  synthesis  should  more  accurately  predict  pilot  vehicle 
performance  than  the  classical  methods  and  this  can  be  shown  by  compar¬ 
ing  root  locus  plots.  Root  locus  plots  for  the  ideal  and  typical  sight 
configurations  are  shown  in  Figures  2  and  3,  respectively,  where  u  shows 
the  unaugmented  system  roots,  and  A-C  indicate  the  root  locations  for 
increasing  levels  of  augmentation.  Notice  that  the  vehicle  short  period 
root  movement  is  completely  different  for  the  two  sight  configurations. 
Since  the  same  vehicle  dynamics  are  present  in  both  cases,  the  only 
difference  between  the  two  cases  would  be  the  higher  order  dynamics 
introduced  by  the  different  sights.  Since  classical  methods  ignore 
these  higher  order  dynamics,  any  augmentation  strategy  based  on  classi¬ 
cal  approaches  may  or  may  not  result  in  the  proper  root  movement.  The 
piloted  vehicle  performance  characteristics  are  different  for  the  two 
sight  configurations,  so  the  higher  order  dynamics  must  obviously  be 
included  in  the  model  in  order  to  accurately  predict  piloted  vehicle 
performance.  The  use  of  the  control  law  synthesis,  therefore,  appears 
justified,  and  a  similar  analysis  will  now  be  done  for  the  more  difficult 
multiaxis  tracking  task. 


Rate,  5(deg/sec)  7.16 


Root  Locus  (Ideal  Sight) 


-12.6  ±  11.  Li 
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2.4  Multiaxis  Tracking  Task 

The  multiaxis  air-to-air  tracking  task  is  very  realistic  of  the 
type  of  task  required  of  fighter  aircraft.  As  is  normally  the  case, 
more  realism  generally  translates  into  more  difficulty,  and  such  is  the 
case  here.  The  difficulty  of  the  multiaxis  tracking  task  arises  because 
of  a  high  degree  of  unsymmetric  coupling  between  the  azimuth  and  eleva¬ 
tion  axes.  The  few  previous  studies  on  multiaxis  tracking  tasks  have 
primarily  been  concerned  with  wings-level,  position  tracking  tasks. ^ 

In  that  type  of  task,  the  dynamics  are  uncoupled,  and  each  axis  can  be 
considered  independently.  The  resulting  models  are  therefore  fairly 
small,  and  very  manageable.  This  section  will  discuss  how  the  pitch- 
only  OCM  is  expanded  to  account  for  the  interaction  between  the  axes  and 
develop  the  parameters  for  the  unique  objective  cost  function  associated 
with  the  multi  axis  task.  Then,  as  in  the  pitch  tracking  task,  the 
augmentation  control  synthesis  will  be  applied  to  determine  what  improve¬ 
ment,  if  any,  is  possible  in  the  piloted  vehicle  performance. 

As  before,  the  system  is  represented  in  state  variable  form  as 
x  =  Ax  +  Bu  and  y(t+0  =  Cx  +  Du.  The  required  matrices  are  obtained 
by  linearizing  the  dynamics  of  an  F-106  aircraft,  using  an  ideal  sight, 
about  a  four-g,  turning  flight  condition,  to  obtain  the  vehicle  pertur¬ 
bations.  The  engagement  parameters  are  presented  in  Table  6,  while  the 
aircraft  state,  control  input,  and  pilot's  observation  vectors  are 
defined  in  Table  7.  A  more  complete  development  of  these  matrices  is 
contained  in  Schmidt's  report,^  and  the  resulting  system  matrices  are 
shown  in  Appendix  A.  A  baseline  optimal  control  model  was  initially 
obtained  as  a  direct  extension  of  the  pitch  tracking  task  model.  The 


Table  6 


Engagement  Parameters 
Target  altitude,  hy  *  10,000  ft. 

Mach  =0.72 

Target  velocity.  Ussy  =  775  ft/sec 

Target  flight  path  angle,  0y  =  0  deg 

Target  bank  angle,  <j>y  =  75.5  deg 

Distance  to  target,  d  =  2000  ft. 

Target/attacker  relative  heading,  au>  =  23.6  deg 

Target  normal  acceleration,  A  =  -128.8  ft/sec 

ZT 

Target  accelerations  (attacker's  coordinates) 

Ay  =  -50.0  ft/sec2 
x 

Ay  =  -  3.0  ft/sec2 

y 

AT  =  -119.0  ft/sec2 
z 

Target  velocities  (attacker's  coordinates) 

Vy  =710  ft/ sec 
x 

Vy  =  78  ft/sec 

y 

Vy  =  -300  ft/sec 
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Table  7 

Linear  Dynamic  Model  Vectors 
State  vector  of  perturbation  variables: 

x  -  £  I.  v »  Vy,  Wy»  d,  »  ^az*  ®*  P»  P»  ] 

A\|»  =  relative  heading  (c-j  -  4^) 

Vj,  Wj  =  target  velocities  (attacker's  coordinates) 
d  =  range  variation 

6EL’  eAZ  =  tar9et  ^ne  sight  angle  components 
0,  o  =  attacker  Euler  angles 
a,  e  =  attacker  angle  of  attack,  sideslip  angle 
p,q,r  =  stability  axis  angular  rates 


Pilot's  observation  vector: 

=  [eEL’  eEL *  eAZ’  eAZ*  XEL’  XEL*  XAZ’  XAZ’  Ay*  ^ 
CEL*  eAZ  =  e^evatlon  and  azimuth  tracking  errors 

xEi»  *  sight  displayed  lead  angle  components 
Ay  =  lateral  acceleration 
4>  =  relative  bank  angle 


Pilot's  control  input  vector  (perturbation  deflections): 
_  T 

Up  _  [<5c  ’  ^ 

p  hst  Ast  Rped 

6r  =  elevator  deflection 
st 

6.  *  aileron  deflection 

Ast 


ped 


=  rudder  deflection 
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parameters  of  the  baseline  model  are  given  In  Table  8.  Notice  that  only 
two  control  Inputs,  the  aileron  and  elevator,  are  Included,  and  the 
significance  of  this  fact  will  be  discussed  later  In  this  section. 
Simulation  data  for  an  F-106  aircraft,  with  several  sight  displays, 
obtained  from  the  large  amplitude,  motion  simulator  at  the  Flight 
Dynamics  Laboratory,  Wrlght-Patterson  Air  Force  Base,  was  used  to  calcu¬ 
late  rms  values  for  various  parameters.  The  rms  performance  character¬ 
istics  of  the  simulation  parameters  were  then  compared  to  the  baseline 
OCM  analytical  results,  and  the  results  are  shown  In  Table  9.  Consider¬ 
ing  that  the  baseline  optimal  control  model  was  essential  a  "single¬ 
axis"  model,  the  values  are  remarkably  close.  The  system's  sensitivi¬ 
ties,  at  least  to  first  order,  were  then  found  by  first  increasing  the 
observation  errors,  and  then  also  increasing  the  neuromuscular  lag  time 
constant,  to  determine  the  changes  in  rms  performance.  The  results  are 
presented  in  Table  10,  and  as  expected,  all  the  rms  values  for  the  states 
and  outputs  increased,  but  it  should  be  noted  that  the  lateral  error 
increased  significantly  more  than  longitudinal  error.  This  tends  to 
Indicate  the  pilot  has  a  greater  difficulty  reducing  lateral  errors  and 
must  therefore  increase  the  attentional  allocation  for  that  axis.  Based 
on  rms  performance  matching,  a  final  model  was  developed. 

The  parameters  for  the  final  OCM  are  given  in  Table  11,  and  there 
are  several  important  points  to  notice.  First,  there  are  three  control 
inputs  instead  of  the  two  control  Inputs  In  the  baseline  model.  The 
rudder  was  added  to  account  for  two  different  tracking  strategies.  The 
first  strategy,  called  the  "rolling"  strategy,  has  the  pilot  first  roll 
the  aircraft  to  null  out  the  azimuth  error,  using  only  ailerons,  and 


— '/•  r* 
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Table  8 
Baseline  OCM 

Observation  vector,  y  *  ^AZ*  ^EL*  ^EL*  ^AZ* 

Objective  function  weights  (both  axes),  Q£  s  16.0,  Q«  *  1.0, 
Observation  thresholds,  *  0.05  deg. 

T*  *  T?  *  0.10  deg. 

Observation  noise  ratio,  -20  dB  for  full  attention 
Fractional  attention,  f^  =  0.5  for  all  observed  variables 
Observation  delay,  t  =  0.2  sec. 

Neuromuscular  lag,  xN  =  0.2  sec.  for  all  control  inputs 

Motor  noise  variance,  varied  to  match  rms  control  inputs 

Control  inputs,  6P  ,6.  only 
Lst  Ast 


XAZ^ 

Q-  -  4.0 


Table  9 

Simulation  vs.  Baseline  Model 
(RMS  Performance) 


eEL(deg) 

'AZ(de9> 

XEL(deg) 

XA2<de9) 

Simulation 

1.09 

0.97 

1.72 

2.58 

Baseline  Model 

1.31 

_ 

0.82 

0.74 

0.34 

VT( ft/sec) 

WT( ft/sec) 

— 

e(deg) 

<t>(deg) 

Simulation 

34 

48.0 

1.09 

4.18 

Baseline  Model 

15.0 

37.0 

0.74 

2.29 

p(deg/sec) 

q( deg/ sec) 

r(deg/sec) 

a(deg) 

Simulation 

4.99 

1.72 

0.92 

1.40 

Baseline  Model 

6.07 

3.04 

1.72 

1.03 

<5F  (in) 
tst 

<$A  (in) 

”st 

6„  (in) 
Rped 

B(deg) 

Simulation 

0.27 

0.22 

0.22 

0.20 

Baseline  Model 

0.28 

0.23 

- 

0.43 

mm'"” 


Table  10 


OCM  Variation  Comparison 
(RMS  Performance) 

Model  A  -  Increased  observation  errors 

Model  B  -  Increased  observation  errors,  plus  Increased  neuromuscular  lag 
constant 


eEL(deg) 

cAZ(de9) 

xEL(deg) 

xA2(deg) 

Baseline  Model 

1.31 

0.82 

0.74 

0.34 

Model  A 

1.60 

1.32 

0.97 

0.52 

Model  B 

1.66 

1.95 

0.34 

0.69 

VT( ft/sec) 

WT( ft/sec) 

e(deg) 

<j.(deg) 

Baseline  Model 

15.0 

37.0 

0.74 

2.29 

Model  A 

22.0 

44.0 

0.97 

3.32 

Model  B 

29.0 

21.0 

1.38 

4.70 

a(deg) 

6(deg) 

HB 

6.  (in) 

Ast 

Baseline  Model 

1.03 

0.43 

0.28 

0.23 

Model  A 

1.09 

0.47 

0.28 

0.27 

Model  B 

1.15 

0.40 

0.24 

0.25 

p(deg/sec) 

q(deg/sec) 

1 

r(deg/sec) 

Baseline  Model 

6.07 

3.04 

1.72 

Model  A 

6.65 

3.09 

1.78 

Model  B 

6.30 

2.58 

1.38 
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Table  11 


Final  Multiaxis  Pilot  Model  Parameters 

Observation  vector,  yT  »  [e^.  PEL*  CAZ*  ^AZ’  XEL*  *El*  *AZ*  *AZ*\’*rel^ 
Objective  function  weightings,  Qe  *  16.0,  Q*  *  1.0,  Qj  *  4.0, 


«  8.0,  QAy  -  .007 


Observation  thresholds,  T£  *  Tx  *  0.05  deg 

T-  =  T*  *  0.18  deg 

C  A 

TAy  =  0.4  ft/sec2 

T,  =  5  deg. 

$ 

Fractional  attention  allocations,  fi  =  0.05  for  e^,  e^L>  x^,  x^L 

fi  =0.15  for  eAZ,  eAZ,  XAZ,  XftZ 
fi  =  0.10  for  Ay,  <j. 

Observation  noise  ratios,  -20dB  for  full  attention 

Observation  delay,  t  =  0.2  sec. 

Neuromuscular  lag,  tn  *  0.33  for  6F  (Q,  *  0.05) 

"  tst  °E 

=  0.23  for  6.  (Q*  =  0.10) 

N  Ast  6A 

tw  =  0.62  for  5D  (Q,  »  0.12) 

N  Rped  6R 


Motor  noise  ratios,  Vn 


wpuTu 


2 


.05(-8dB)  for  5F 
bst 


p„  *  .02{-12dB)  for  6. 
u  Ast 

p,  *  .05(-8dB)  for  6„ 
u  Rped 
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then  correct  the  elevation  errors,  using  the  elevator,  with  small 
azimuth  corrections,  using  the  ailerons.  A  second,  or  "pointing", 
strategy  could  also  be  used,  however,  whereby  the  pilot  uses  ailerons, 
elevator,  and  rudder  to  continually  "point"  at  the  target.  This  latter 
strategy  requires  a  considerable  amount  of  "cross  controlling"  of  the 
control  surfaces,  however.  The  simulation  results,  shown  in  Table  9, 
indicate  that  sideslip,  6,  is  fairly  small,  and  hence  a  "rolling"  stra¬ 
tegy  was  probably  employed.  From  my  own  personal  experiences  as  an  F-4 
fighter  pilot,  I  tend  to  agree  that  the  majority  of  tracking  is  done 
using  the  "rolling"  strategy,  but  that  as  the  target  range  gets  smaller, 
or  target  motion  increases,  a  "pointing"  strategy  must  often  be  used  to 
maintain  the  sight  on  the  target.  Therefore,  the  three  inputs  were 
included  to  allow  for  all  possibilities.  The  second  point  worthy  of 
notice  is  the  change  in  the  pilot's  observation  vector,  where  the 
lateral  acceleration,  aT,  and  relative  bank  angle,  <j>y,  have  been  added. 
Lateral  acceleration  was  included  because  test  pilots  indicated  that 
pilots  tend  to  minimize  lateral  acceleration,  even  if  larger  errors 
result.  Again,  from  personal  experience,  I  can  verify  this  fact.  The 
inclusion  of  the  relative  bank  angle  seems  appropriate,  since  fighter 
pilots  use  target  relative  bank  angle  to  help  predict  the  target's 
movement.  While  the  short  range  tracking  solutions  simulated  in  this 
research  tend  to  keep  this  angle  small,  it  could  be  very  important  for 
future  studies.  The  analytical  performance  predictions  were  then 
compared  to  simulation  results,  and  both  are  shown  in  Table  12.  The 
double  asterisks  indicate  that  the  final  model  provided  the  best  cor¬ 
relation  between  analytical  and  simulation  results,  while  the  single 
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cEL(deg) 


XAZ(deg) 


Simulation 

1.09 

0.97 

1.72 

2.58 

Analytical 

1.78 

1.72 

1 .38** 

0.63* 

VT(ft/sec) 

Wj(ft/sec) 

©(deg) 

<t>(deg) 

Simulation 

34.0 

48.0 

1.09 

4.18 

Analytical 

25.0* 

63.0 

0.97* 

3.55* 

Simulation 

Analytical 


<sF  (in) 

<5a  (in) 

Lst 

Ast 

0.27 

0.22 

0.23 

0.15 

o  (deg) 


p(deg/sec) 

q(deg/sec) 

- - 

n( deg/sec) 

B(deg) 

Simulation 

Analytical 

4.99 

5.21** 

1.72 

2.58** 

_ 

0.92 

0.57** 

0.20 

0.15** 

— 

Tj^(sec) 

xN  (sec) 

na 

Vsec) 

- - 

Ay(ft/sec  ) 

Simulation 

- 

- 

- 

0.67 

Analytical 

.33 

.28 

_f _ 1 

0.78** 

l  .  «u»i 


err*  - 
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asterisk  means  one  other  model  was  better.  Since  only  one  pilot  flew 
the  simulations,  and  only  one  simulation  run  met  the  required  engagement 
parameters  and  sight  configuration,  further  correlation  would  be  point¬ 
less.  Therefore,  the  final  optimal  control  model  will  be  used  as  the 
system  model  for  applying  the  augmentation  control  synthesis. 

The  augmented  control  synthesis  for  the  multi  axis  tracking  task 
follows  a  similar  procedure  to  the  pitch  tracking  task.  Because  the 
"rolling"  strategy  was  assumed  to  be  used,  no  control  gains,  however, 
were  put  on  the  rudder.  The  analytical  predictions  of  rms  performance, 
using  the  augmented  system,  are  given  in  Table  13,  and  a  significant 
reduction  in  tracking  errors,  particularly  in  azimuth  errors,  is 
indicated.  This  indicates  that  the  augmented  system  makes  the  largest 
impact  in  the  areas  where  the  pilot  has  the  greatest  difficulty.  It 
should  also  be  noted  that  the  stick  and  rudder  deflections  have  been 
reduced,  which  indicates  a  reduction  in  pilot  workload.  Therefore,  the 
augmentation  control  synthesis  yields  a  marked  improvement,  at  least 
analytically,  in  the  piloted  vehicle  performance.  Some  small  changes 
will  now  be  added  to  the  optimal  control  model  to  more  realistically 
simulate  the  multiaxis  air-to-air  tracking  task. 

2.5  Model  Modification 

The  primary  objective  of  this  research  is  to  simulate  a  multiaxis 
air-to-air  tracking  task  and  use  the  simulation  data  to  validate  the 
augmentation  control  synthesis,  and  if  necessary,  to  change  the  OCM 
parameters  to  better  predict  the  piloted  vehicle  performance.  In  an 
attempt  to  more  realistically  simulate  the  task,  two  changes  have  been 
made  to  the  final  OCM.  First,  the  target  range  is  set  to  a  constant 
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Unaugmented 

Augmented 


Table  13 

Augmented  System  Performance 
(Analytical  RMS  Results) 


eEL(deg) 

eAZ(deg) 

xEL(deg) 

*AZ(deg) 

1.78 

1.72 

1.38 

0.63 

0.54 

0.20 

0.33 

0.10 

VT( ft/sec) 

W^.(  ft/sec) 

o(deg) 

4-  ( deg) 

Unaugmented 

25.0 

63.0 

0.97 

3.55 

Augmented 

4.0 

15.0 

0.17 

— 

0.48 

Augmented 


6p  (in) 

6*  (in) 

hst 

st 

0.23 

0.15 

0.19 

1 

0.11 

p( deg/sec) 

q(deg/sec) 

n(deg/sec) 

6(deg) 

Unaugmented 

5.21 

2.58 

0.57 

0.15 

Augmented 

1.20 

0.97 

0.21 

0.06 

!  T^(sec) 

tn  (sec) 

WA 

(sec) 

nR 

Ay( ft/sec2) 

Unaugmented 

0.33 

0.28 

0.62 

0.78 

Augmented 

0.33 

0.28 

0.62 

0.25 
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2000  feet.  This  was  done  primarily  to  make  the  analysis  of  data  easier. 
The  second  change  introduces  perturbations  to  the  constant  bank,  four-g 
turn  used  to  initially  develop  the  target  motion.  Target  perturbations 
are  introduced  by  assuming  the  pilot  will  maintain  a  constant  four-g 
turn,  but  his  bank  angle  will  perturbate  about  the  steady-state  value 
of  75.5°.  To  determine  how  these  bank  angle  changes  will  affect  the 
target  motion,  the  assumptions  used  to  develop  the  target's  motion  must 
be  studied.  The  values  of  Vy  and  Wy,  in  the  x-vector,  are  the  target’s 
velocity  perturbations,  in  the  attacker's  coordinate  system.  The  equa¬ 
tions  governing  them  assumed  vT  and  wT  ,  the  target's  velocity  pertur- 

tT  't 

bations  from  the  steady-state  values,  in  the  target's  coordinate  system, 

are  zero.  When  the  bank  angles  vary  about  the  steady-state  target  bank 

angle,  this  assumption  is  no  longer  true,  and  vT  and  wT  must  be  added 

T  .  T  . 

to  the  "vT  and  wx  equations.  However,  the  values  of  vT  and  wT  must  be 
T  T  tT  *T 

transformed  from  the  target's  coordinate  system  to  the  attacker's  coor¬ 
dinate  system  before  they  can  be  added  together.  The  coordinate  trans¬ 
formation  to  accomplish  this  is  fairly  straight  forward,  and  a  more 
detailed  development  is  shown  in  Appendix  B.  Now,  it  is  necessary  to 


show  how  to  solve  for  vT  and  wT  . 

T  T 

The  target  motion  perturbation  equations  are  given  by  Vy  = 

(g  cos  <i>i ) 4>i-  and  Wy  =  (-g  sin  4>j)#y.  wi.ere  ^  =  75.5°,  the  steady- 
state  target  bank  angle,  and  =  the  target  perturbation  bank  angle. 
Therefore,  the  problem  is  to  solve  for  a  perturbation  bank  angle  that 


realistically  simulates  a  real  target's  motion.  First,  $y  was  assumed 

to  be  small,  so  that  qT  and  rT  ,  the  target  perturbation  Euler  angles, 

tT  tT 

are  assumed  equal  to  zero.  Therefore,  the  target  perturbations 


(1  93^ 

equations  are  given  by  ij<y  =  0,  Qy  -  -.161  $y,  and  +  .161  Q-p 

The  value  of  pT  is  then  calculated  using  pT  =  -  —  pT  + 

iT  1  tr  1  A 

where  tr  is  the  target  vehicle  roll  time  constant,  and  6^  is  the  aileron 

control  surface  deflection,  6.,  is  calculated  by  using  a  sum  of  sinusoids 
13  * 

using  6.(t)  =  ^  A,  sin(w.t+<t>. ) ,  with  the  variables  defined  and  tabu- 

a  i=1  11 

lated  in  Appendix  c.  Once  the  value  of  6^  is  computed,  the  target  per¬ 
turbation  equations  for  py,  ©y,  and  4>j  are  integrated  to  determine  ij>y, 

which  in  turn  is  used  to  compute  vT  and  wT  .  These  values  are  then 

T  T 

transformed  into  the  attacker's  coordinate  system  and  added  to  the  \fy 
and  Wy  equations,  respectively.  One  last  note;  the  value  of  xR,  the 
target's  vehicle  roll  constant,  was  set  at  2.0,  which  is  the  approximate 
T-38  roll  time  constant  at  the  given  flight  conditions.  These  small 
changes  to  the  model  result  in  a  much  more  realistic  tracking  problem. 


2.6  Summary 

The  optimal  control  model  (OCM)  is  found  to  be  an  extremely  valuable 
tool  for  predicting  piloted  vehicle  performance  when  higher  order  vehicle 
dynamics  are  present.  The  OCM,  once  the  objective  cost  function  for  the 
particular  task  is  discovered,  is  a  fairly  simple  and  straight  forward 
method  of  analyzing  piloted  vehicle  performance.  The  comparison  between 
the  analytical  and  simulation  results  provides  excellent  correlation. 

A  systematic  augmentation  control  synthesis  procedure  for  computing  the 
augmentation  gains  is  then  applied  to  the  OCM  to  improve  the  piloted 
vehicle  performance.  The  control  synthesis  has  been  applied  to  a  pitch 
tracking  task,  using  aircraft  with  two  sight  configurations,  whose 
dynamics  result  in  higher  order  system  dynamics.  The  analytical  predic¬ 
tions  of  improved  piloted  vehicle  performance,  when  augmentation  is 
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added,  must  still  be  verified  experimentally,  but  a  vast  improvement  in 
performance  is  theoretically  possible.  A  similar  analysis,  using  the 
more  realistic  multiaxis  tracking  task,  shows  similar  results,  but  it, 
too,  requires  experimental  verification,  which  is  the  primary  objective 
of  this  research.  To  provide  for  a  realistic  simulation  task,  an 
enhanced  degree  of  target  motion  has  been  added  to  the  OCM.  The  next 
chapter  will  provide  a  detailed  analysis  of  how  the  multiaxis  tracking 
task  will  be  simulated. 


CHAPTER  3.  SIMULATION  SYSTEM 


3.1  Introduction 

The  Purdue  University  Aeronautical  and  Astronautical  Engineering 
(A&AE)  Flight  Simulation  Laboratory  provides  the  basis  for  simu¬ 
lating  the  multi-axis  air-to-air  tracking  task.  This  fixed  base  simula- 
tion  was  originally  developed  in  1978,'  ;  and  has  been  used  for  several 
research  simulations  here  at  Purdue  University.  A  more  complete 
discussion  of  the  system  can  be  found  in  Reference  8,  but  this  chapter 
will  first  present  a  brief  overview  of  the  hardware  components  used  in 
the  simulation,  emphasizing  their  functions  in  the  simulation,  not  their 
specifications.  Then,  the  development  of  the  software  to  simulate  the 
task  will  be  detailed. 

3.2  System  Hardware 

The  hardware  used  in  the  simulation  can  be  divided  into  four  main 
areas:  the  computer,  the  display  mechanisms,  the  control  input  apparatus, 
and  the  input-output  (I/O)  devices.  This  section  will  discuss  the 
function  and  importance  of  each  of  these  areas. 

The  heart  of  any  computer  simulation  is  the  computer,  since  it  not 
only  must  calculate  the  various  parameters  used  to  update  the  display, 
but  also  acts  as  the  command  center  for  controlling  the  other  pieces  of 
hardware.  The  General  Automation  SPC-T5/45  minicomputer  performs 
these  functions  and  has  full  I/O  capability  with  the  disk,  plus 
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a  complete  range  of  real  time  clock  functions,  which  makes  the  mini¬ 
computer  well  suited  for  this  type  of  simulation.  The  minicomputer  is 
also  capable  of  performing  I/O  operations  with  several  other  I/O  devices, 
which  will  be  detailed  later  in  this  section.  The  minicomputer  obviously 
calculates  various  parameters,  which  will  be  fully  explained  in  Section 
3.3,  but  it  also  acts  as  the  workhorse  for  updating  the  graphics  display. 
The  GA  SPC  15/45  uses  a  display  listing,  which  is  a  defined  storage 
location,  where  various  display  "entities"  are  stored.  The  minicomputer 
then  executes  various  graphics  instructions  that  individually  change  or 
move  the  various  entities,  and  these  binary  plotting  instructions  are 
then  sent  to  the  display  mechanisms.  The  six  entities  and  set  of  possi¬ 
ble  graphics  instructions  used  in  the  simulation  are  presented  in 
Appendix  E. 

The  display  mechanisms  consist  of  two  pieces  of  equipment:  the 
QVEC  2150  vector  generator  and  the  Hewlett-Packard  1310A  cathode  ray 
tube  (CRT).  The  vector  generator  is  the  display  processor,  which 
accepts  the  incoming  binary  plottings  instructions  and  calculates  the 
positioning  voltages  and  on/off  signals  used  to  update  the  CRT  display. 
The  vector  generator  also  has  direct  access  to  the  display  listing  con¬ 
tents  of  the  minicomputer,  so  that  the  computer  execution  time  for 

updating  the  CRT  electron  beam  instructions,  which  must  be  done  about 

(8) 

30  times  per  second  to  prevent  flickering,  is  kept  to  a  minimum.  ' 

The  Hewlett-Packard  CRT  is  an  ll"xl5"  screen  with  a  50  lines/inch 
resolution,  and  the  execution  time  for  drawing  the  beams  is  essentially 
zero. ^  These  two  pieces  of  hardware  combined  to  form  an  effective 


method  for  displaying  a  simulation  display. 
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Another  important  piece  of  simulation  hardware  is  the  apparatus 
for  obtaining  the  control  inputs.  The  apparatus  includes  a  nonmoving 


stick  to  obtain  elevator  and  aileron  control  inputs  and  a  nonmoving  bar 
to  obtain  rudder  inputs,  and  is  similar  to  some  newer  fighter  aircraft 
designs,  such  as  the  F-16.  When  a  force  is  applied  to  the  stick  or 
rudder,  resistance  strain  gages  measure  a  voltage  that  is  proportional 
to  the  applied  force. ^  The  maximum  output  voltage  is  +  10  v, 
but  the  amount  of  force  required  to  obtain  the  maximum  voltage  can  be 
changed  by  varying  the  control  gains,  which  is  done  in  the  software, 
and  will  be  further  explained  in  Section  3.3.  Trim  knobs  are  also 
included  that  can  adjust  the  zero-force  voltage  by  +1  volt.  The  method 
of  reading  these  voltages,  and  then  converting  them  to  a  control  surface 
deflection,  is  again  accomplished  by  the  software,  and  will  be  explained 
in  Section  3.3. 

The  final  hardware  component  area  is  input -output  devices,  and  in 
addition  to  the  disk  unit  in  the  minicomputer,  includes  a  card  reader, 
a  Superbee  teletype  terminal,  a  line  printer,  and  the  VERSATEC  electro¬ 
static  plotter/printer.  The  minicomputer  disk  is  primarily  used  to 
store  working  computer  programs  and  act  as  a  simulation  data  storage 
unit.  The  data  storage  on  the  disk  provides  an  easy  method  of  catalog¬ 
ing  the  simulation  data  and  still  allows  an  easy  method  of  transferring 
the  data  to  more  permanent  storage  devices,  such  as  magnetic  tape.  The 
elements  of  various  matrices  used  to  update  the  simulation  parameters 
can  be  inputted  into  the  computer  program  by  either  the  card  reader  or 

by  using  a  storage  location  on  the  minicomputer  disk.  The  latter  method 
was  chosen  for  this  simulation  because  the  elements  of  these  matrices 
are  constant.  If  the  elements  varied  between  runs,  reading  this 


information  on  the  card  reader  would  be  more  convenient.  The  teletype 
terminal  could  also  be  used  to  input  information,  but  the  primary  use 
is  for  editing  programs.  The  teletype  terminals  can  also  be  useful 
for  rapidly  looking  at  the  simulation  output,  while  the  line  printer 
provides  a  more  permanent  printed  record.  The  electrostatic  printer/ 
plotter  is  especially  useful  because  the  simulation  data  can  be 
plotted  almost  immediately.  It  can  also  be  used  to  pring  the  same 
information  as  the  line  printer.  A  more  complete  description  of  how 
each  of  these  I/O  devices  are  used  for  this  experiment  is  presented  in 
Chapter  4. 

The  hardware  used  for  the  multi axis  tracking  task  simulation 
provides  all  the  equipment  to  successfully  accomplish  the  simulation. 

The  GA  minicomputer,  combined  with  the  vector  generator  and  cathode 
ray  tube,  calculates  the  appropriate  simulation  parameters  and  then 
updates  the  display.  The  various  input/output  devices  allow  for  a 
wide  range  of  inputting  information,  storing  and  transferring  data, 
and  presenting  results.  While  the  hardware  is  the  backbone,  or  frame¬ 
work,  of  the  simulation,  the  software  is  the  "brains"  of  the  simulation 
and  will  be  discussed  next. 

3.3  System  Software  Development 

The  computer  software  which  accomplishes  the  multi  axis  tracking 
task  simulation  can  be  divided  into  four  phases:  the  initialization 
and  precomputation  phase,  the  vector  update  phase,  the  display  update 
phase,  and  the  data  storage  and  transfer  phase.  A  computer  flow  chart 
is  shown  in  Figure  3  and  provides  a  overview  of  the  overall  programming 
sequence.  The  actual  simulation  computer  listing  is  given  in  Appendix  F 
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and  also  includes  variable  definitions  and  subroutine  explanations. 

This  section  will  outline  the  development  of  each  phase  of  the  software 
and  give  an  explanation  of  why  certain  methods  were  chosen. 

The  initialization  and  precomputation  phase  begins  by  setting  up 
a  data  file  to  allow  inputting  of  information  through  the  teletype 
terminals.  The  next  step  is  to  "read"  in  the  values  of  the  appropriate 
matrices.  Once  this  is  accomplished,  the  display  area  is  initialized, 
which  means  that  the  name  and  dimension  of  the  storage  area  containing 
the  display  list  is  determined.  Also,  two  variables,  IAUG  and  DT,  are 
input  via  the  Superbee  teletype  terminal.  The  variable,  IAUG,  is  a 
logical  variable  that  determines  if  the  augmentation  system  is  included 
in  the  simulation  run  or  not.  If  IAUG  =  0,  no  augmentation  is  included, 
while  any  other  value  will  provide  system  augmentation.  The  parameter, 
DT,  is  the  iteration  step  size,  At,  and  its  value  and  importance  will 
be  explained  later  in  this  section.  These  variables  are  input  in  this 
manner  to  allow  for  changing  the  parameters  once  the  program  is  com¬ 
piled.  The  value  chosen  for  DT  will  be  further  explained  later  in  this 

section.  The  next  step  in  this  phase  is  to  calculate  the  augmentation 

•  AAt 

matrices,  AaUg  and  Baug,  and  the  state  transition  matrices,  e  and 

(  e  dt)B,  where  A  and  B  represent  either  the  augmented  or  unaugmented 
Jo 

system  matrices.  In  the  actual  programing,  only  Agug  and  Baug,  defined 
as  AAUG  and  BAUG,  are  used,  because  when  IAUG  *  0,  AAUG  and  BAUG  are 
set  equal  to  the  A  and  B  matrices,  respectively,  while  when  IAUG  ^  0, 
AAUG  and  BAUG  are  calculated  using  the  control  synthesis  relationships 
of  A =  A-Bkv  and  Bail_  *  B-Bk„.  The  calculation  of  the  state  transi- 

ouy  a  aug  u 

tion  matrices  is  done  by  a  series  approximation  and  the  details  are 


presented  in  Appendix  D.  Finally,  prior  to  beginning  the  iteration 
cycle,  all  the  initial  values  are  set. 

The  most  important  phase  of  a  computer  simulation  is  the  vector 
update  phase,  where  the  state  and  observation  vectors  are  updated  at 
each  iteration  step,  At.  Before  detailing  the  method  used,  a  brief 
discussion  of  why  the  state  transition  method  was  chosen  will  be  given. 
The  two  important  parameters  for  updating  the  state  and  observation 
vectors  are  speed  and  accuracy.  The  speed  is  important,  because  the 
display  is  only  updated  once  during  each  iteration  cycle,  and  a  total 
cycle  time  greater  than  about  .1  second  will  not  adequately  simulate 
the  task.  This  is  because  the  minicomputer  obtains  discrete  solutions, 
which  must  reasonably  represent  the  aircraft's  highest  frequency  mode, 
and  a  value  of  At  less  than  .1  second  provides  a  reasonable  speed  of 
display  updating. The  accuracy  is  obviously  important,  because 
the  values  of  the  state  and  observation  vectors  are  used  to  update  the 
display.  The  original  A&AE  minicomputer  system  was  designed  to  simulate 
aircraft  motion  using  six  degree  of  freedom,  nonlinear  equations  of 
motion.  A  fourth  order  Runga-Kutta  integration  routine  was  found  to 
have  both  the  required  speed  and  accuracy  to  meet  the  original  simula¬ 
tion  needs. ^  Several  other  simulations  also  used  this  method, 
even  though  linear  equations  of  motion  were  used.  Therefore,  it  seemed 
reasonable  to  incorporate  the  fourth  order  Runga-Kutta  routine  for  the 
multiaxis  tracking  task.  A  problem  arose,  however,  because  the  itera¬ 
tion  cycle  time.  At,  was  nearly  0.2  seconds,  which  was  clearly  unaccept¬ 
able.  The  increase  in  cycle  time  resulted  from  the  added  complexity  of 
the  task,  since  the  previous  simulations  were  primarily  pitch  tracking 
tasks,  and  the  state  and  observation  vector  dimensions  were  less  than 
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half  the  dimension  of  the  multiaxis  case.  Therefore,  the  state  transi¬ 
tion  method  was  chosen  as  a  second  candidate.  The  state  transition 
method  calculates  the  new  state  values  using 

x(t+At)  =  e AAt  x(t)  +  ([  eAxdT)  Bu(t)  .  (3.3.1) 

J  o 

This  method  will  calculate  exact  answers  if  the  inputs,  u[t),  are  held 
constant  over  the  interval,  At.  Since  the  inputs  are  only  sampled  once 
during  each  cycle,  they  are  not  constant  over  the  interval.  However, 
the  fourth  order  Runga-Kutta  method  uses  the  inputs  in  a  similar  method, 
so  accuracy  of  either  method  seems  reasonable.  The  state  transition 

A. f  fAt  A 

matrices,  e^1  and  (  e  dt)B,  are  constant,  once  the  interval  time, 

•*0 

At,  is  specified,  so  they  can  be  calculated  prior  to  the  iteration  cycle. 
The  savings  in  iteration  computation  time  was  considerable,  and  the 
total  iteration  cycle  time  was  reduced  to  approximately  0.076  seconds. 

One  drawback  is  that  the  cycle  time.  At,  cannot  vary,  and  as  will  be 
seen  later,  the  data  transfer  method  requires  an  additional  0.0125 
seconds,  but  only  every  seventh  iteration.  Therefore,  the  At  for  each 
iteration  must  account  for  this  added  cycle  time.  Since  the  overall  At 
needs  to  be  greater  than  0.0887,  the  At  still  remains  within  acceptable 
limits,  so  no  serious  problem  results.  The  state  transition  method, 
then,  was  found  to  have  both  the  required  speed  and  accuracy,  and  is 
used  to  update  the  state  vector. 

The  vector  update  phase  uses  the  control  input  values  to  calculate 
the  new  values  of  the  state  and  observation  vectors,  and  then  calculates 
the  target  motion  by  integrating  the  target  motion  perturbation  equations. 
The  control  input  apparatus  was  explained  in  Section  3.2,  and  a  FORTRAN 
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subroutine,  DASIN,  converts  the  analog  voltages,  which  are  proportional 

to  the  applied  stick  forces,  into  a  digital  voltage  by  averaging  four 

voltage  readings  over  a  1  msec  time  period.  This  voltage  is  then 

multiplied  by  a  control  device  gain  (in/volt)  that  converts  the  stick 

and  rudder  pedal  forces  into  control  surface  deflections  (in).  The 

gains  were  obtained  experimentally  to  yield  a  reasonable  "feel"  of  the 

aircraft  and  set  at  1.0,  0.5,  and  0.5  for  elevator,  aileron,  and  rudder. 

These  control  input  deflection  values  make  up  the  u  vector  and  are  used 

P 

to  calculate  the  new  values  for  the  state  and  observation  vectors,  using 

x(t+At)  =  e^  x(t)  +  (f  e  dx)Bu(t)  (3.3.2) 

■'o 

y(t+At)  =  Cx(t)  +  Du(t)  .  (3.3.3) 

The  next  step  in  the  vector  update  phase  is  to  calculate  the  target 
motion.  The  perturbation  equations  governing  the  target  motion  were 
developed  in  Section  2.5  and  are  repeated  here. 


6A 

(3.3.4) 

5 

1 

II 

(3.3.5) 

*T  =  PT  +  -161 

°T  • 

(3.3.6) 

The  perturbation  target  aileron  deflections,  6.,  were  modeled  using  a 

13 

sum  of  sinusoids,  where  6.(t)  =  T  A.  sin  (w.t+<j>.),  and  the  values  of 

a  i=l  1  1  1 

Aj,  Wj,  and  ^  are  given  in  Appendix  C,  and  were  chosen  to  simulate  some 
realistic  target  motion.  Once  <5^(t)  is  calculated,  the  perturbation 
equations  are  integrated  using  the  fourth  order  Runga-Kutta  integration 


routine,  to  obtain  4>j,  which  is  used  to  calculate  vT  and  wT  .  As 

T  T 

explained  in  Section  2.5,  vT  and  wT  must  be  transformed  from  the 

‘t  't 

target's  coordinate  system  into  the  attacker's  coordinate  system.  The 
actual  transformation  matrix  development  is  shown  in  Appendix  D,  and 
the  actual  coding  in  the  computer  program  was  chosen  to  minimize  compu¬ 
tation  time.  Once  the  state  and  observation  vectors  and  the  target 
motion  have  been  calculated,  the  simulation  display  can  be  updated. 

The  next  phase  of  the  simulation  program,  where  the  simulation 
display  is  updated,  plays  an  important  part  in  the  simulation,  because 
the  pilot's  strategy,  or  selection  of  control  inputs,  is  based  on 
reacting  to  the  visual  display.  It  is  especially  important  for  this 
system  because  the  pilots  only  feedback  of  his  control  input  is  through 
the  display.  The  simulation  displays,  as  depicted  in  Figure  4,  consists 
of  a  stationary  reference  cross,  a  sight  symbol,  a  target  symbol,  and 
three  horizon  lines.  As  mentioned  in  Section  3.2,  the  minicomputer 
stores  these  entities  in  the  display  list,  and  various  FORTRAN  graphics 
instructions  place  them  in  the  proper  location  on  the  CRT.  The  list 
of  entities  and  the  set  of  possible  graphics  instructions  are  contained 
in  Appendix  E.  The  reference  cross,  which  is  set  when  the  display  area 
is  initialized,  represents  the  attacking  aircraft's  forward  velocity 
vector,  and  as  mentioned,  remains  fixed  on  the  screen.  The  other  five 
entities  are  moved  and  rotated  on  the  CRT  using  the  graphics  instructions, 
based  on  the  state  and  output  vectors  and  target  motion  values.  For 
example,  the  sight  and  target  locations  are  positioned  on  the  CRT  using 
the  instruction,  GVECT,  and  the  values  and  g^,  which 

are  the  Y(5),  Y{7),  X(4),  and  X(5)  values,  respectively.  In  addition, 
the  target  symbol  is  rotated,  based  on  4>j,  which  was  found  in  the  target 
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motion  calculations.  The  values  of  A  and  B  are  angles,  measured  in 
radians,  that  must  be  converted  into  raster  units  by  multiplying  them 
by  scale  factors,  to  determine  the  actual  "distance"  on  the  CRT  the 
symbols  are  located.  These  scaling  factors  are  somewhat  arbitrary,  and 
were  chosen  as  300.0,  based  on  experimentation,  to  be  able  ic 
provide  "reasonable"  movement  of  the  sight  and  target.  My  personal 
experiences  as  an  F-14  fighter  pilot,  with  included  many  hours  of 
air-to-air  combat  flying,  were  used  to  make  the  display  as  realistic 
as  possible.  The  horizon  lines  are  included  to  give  the  pilot  an  easy 
visual  cue  of  his  bank  angle,  which  is  calculated  by  adding  the 
attacker's  perturbation  bank  angle,  X(7),  to  75.5°,  the  steady-state 
bank  angle.  The  lines  are  drawn  using  trigonometry,  and  the  appropriate 
graphics  instructions,  to  position  three  equidistant  horizon  lines. 

Once  the  simulation  display  has  been  updated,  the  data  can  be  stored 
and  transferred. 

The  final  phase  of  the  computer  software,  which  is  the  data  storage 
and  transfer  phase,  is  important  not  only  for  future  evaluation  of  the 
simulation,  but  also  because  the  method  used  could  possibly  interfere 
with  the  iteration  cycle  time.  The  actual  variables  that  are  stored  at 
the  end  of  each  iteration  cycle  are  the  time,  T,  the  twelve  elements  in 
the  state  vector,  x(t),  the  first  nine  elements  of  the  observation 
vector,  y(t),  the  three  control  inputs,  up,  the  target's  aileron  deflec¬ 
tion,  <5^,  and  the  target  perturbation  bank  angle,  These  values  are 
initially  stored  in  a  28  by  7  storage  matrix,  XBUF  at  the  end  of  each 
iteration.  At  the  end  of  every  seventh  iteration,  a  FORTRAN  subroutine, 
called  FSDIO,  moves  the  entire  storage  matrix  to  a  buffer  location. 
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where  the  disk  can  store  the  data  while  the  minicomputer  simultaneously 
continues  the  simulation  computations.  The  movement  of  data  from  the 
storage  location  to  the  buffer  location  requires  approximately  0.0125 
seconds,  which  was  mentioned  earlier  in  this  section.  This  method  was 
chosen  because  it  is  the  fastest  method  of  data  transfer  available  on 
this  system,  and  because  it  interferes  the  least  with  the  simulation 
computations,  which  keep  iteration  cycle  time  low.  The  28x7  matrix  is 
used  because  FSD10  is  capable  of  transferring  400  words  of  memory,  or 
200  values,  and  smaller  amounts  of  data  transfer  do  not  appreciably 
reduce  transfer  time.  At  the  end  of  the  simulation  run,  the  data  is 
moved  to  a  labelled  storage  area,  where  it  is  available  for  data  analy¬ 
sis.  Also  at  this  time,  any  additional  values  can  be  postcomputed  and 
added  to  the  permanent  storage  location.  For  this  simulation,  the 
equivalent  augmented  control  inputs  are  postcalculated.  More  informa¬ 
tion  about  the  cataloguing  and  analysis  of  data  can  be  found  in  Section 
4.4. 


2.4  Summary 

The  simulation  system  at  Purdue  University  used  to  simulate 
the  multiaxis  tracking  task  provides  a  realistic  simulation  platform. 
The  GA  SPC- 15/45  minicomputer  is  the  heart  and  workhorse  of  the  system 
and  not  only  calculates  the  appropriate  parameters,  but  also  controls 
the  display  mechanisms  and  input/output  devices.  The  nonmoving  control 
input  apparatus  uses  strain  gage  measurements  to  convert  control  stick 
and  rudder  pedal  forces  into  control  surface  deflections.  The  software 
development  strived  for  both  accuracy  and  speed,  in  order  to  make  the 
simulation  as  realistic  as  possible.  The  state  transition  method  was 
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chosen  to  update  the  state  vector  for  just  these  reasons.  The  simula¬ 
tion  display  updates  five  display  entities,  based  on  the  calculated 
values  in  the  state  and  observation  vectors  and  the  target  motion 
calculations.  The  combination  of  the  hardware  and  software  provide  for 
a  highly  realistic  simulation  of  the  multiaxis  air-to-air  tracking  task. 
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CHAPTER  4.  EVALUATION  PROCESS  AND  RESULTS 

4.1  Introduction 

The  purpose  of  any  research  evaluation  is  to  use  some  method  of 
determining  the  validity  of  the  research  hypotheses.  The  primary 
hypothesis  for  this  research  is  that  the  augmentation  control  synthesis 
will  result  in  improved  piloted  vehicle  performance,  while  a  secondary 
hypothesis  is  that  RMS  statistics  obtained  from  a  simulation  can  be 
compared  to  RMS  statistics  for  an  analytical  study  to  determine  if  any 
OCM  parameters  should  be  changed.  The  method  for  testing  both  hypothes¬ 
es  involves  simulating  a  multiaxis  tracking  task  for  both  the  unaugment¬ 
ed  and  augmented  configurations.  Several  pilots  take  part  in  the 
simulation  so  that  average  RMS  statistics  can  be  calculated  for  various 
parameters.  In  addition,  the  pilot's  subjective  ratings  of  the  augmen¬ 
tation  system  will  help  determine  the  acceptability  of  the  augmented 
system.  While  every  researcher  hopes  the  evaluation  will  substantiate 
the  hypotheses,  the  process  will  also  help  determine  the  validity  of 
any  assumptions  and  hopefully  give  some  clues  towards  ways  of  improving 
the  model.  This  chapter  will  first  examine  some  of  the  model  assumptions 
and  then  explain  the  evaluation  process.  An  outline  of  the  simulation 
and  RMS  calculation  procedures  will  then  be  given.  Finally,  a  summary 
of  simulation  and  analytical  results  will  be  tabulated  and  discussed. 
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4.2  Assumptions 

This  section  will  examine  two  model  assumptions  and  explain  how 
they  were  handled  in  the  simulation.  The  initial  assumption  says  the 
target  perturbation  bank  angle,  which  introduces  the  target  pertur¬ 
bation  motion,  is  a  small  angle.  This  assumption  was  made  to  simplify 
the  transformation  matrix  computations.  The  calculation  of  <f»y  is  out¬ 
lined  in  Section  3.3,  and  the  parameters  and  methodology  was  chosen  to 
obtain  a  mean  and  standard  deviation  of  approximately  0  and  5  degrees, 
respectively,  for  a  90  second  simulation  run.  The  actual  values  for 
each  individual  simulation  run  can  be  found  in  Appendix  H  and  the 
average  simulation  values  were  -.38  and  5.25.  Therefore,  the  small 
angle  assumption  for  <t>y  is  assumed  to  be  valid  for  the  simulation 
evaluation. 

The  second,  and  more  important,  assumption  requires  a  well  trained, 
highly  motivated  individual  to  be  accomplishing  the  task.  The  parame¬ 
ters  of  the  OCM  and  the  control  synthesis  are  both  based  on  this  premise, 
so  any  comparison  between  analytical  and  simulation  results  must  insure 
this  assumption  is  true.  The  initial  problem  was  to  find  individuals 
who  could,  after  some  minimum  amount  of  training,  adequately  track  the 
target  to  meet  the  well  trained  assumption.  The  criteria  for  determin¬ 
ing  a  well  trained  pilot  was  the  ability  to  keep  the  sight  symbol 
touching  some  part  of  the  target  symbol  for  at  least  90  seconds.  If  an 
individual  could  not  attain  this  proficiency  within  about  a  30  minute 
initial  training  period,  they  were  not  used  for  the  simulation  evalua¬ 
tion.  Several  individuals  tried  the  task,  but  only  four  people  were 
able  to  meet  the  criteria.  Three  of  the  individuals  are  U.S.  Air  Force 
pilots,  and  additional  information  about  all  four  individuals  can  be 
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found  in  Appendix  H.  After  this  initial  training  session,  several  30 
second  data  runs  were  made  to  not  only  improve  proficiency,  but  to  help 
determine  how  long  each  simulation  run  should  be.  More  information 
about  this  will  be  presented  in  Section  4.3.  By  the  time  the  data  was 
collected,  each  individual  could  be  considered  well  trained  in  the  task. 
There  is  no  easy  objective  criteria  for  determining  a  highly  motivated 
individual,  but  each  subject  is  a  graduate  student  who  works  in  the 
controls  area  of  the  Department  of  Aeronautics  and  Astronautics  at 
Purdue  University,  and  therefore  understands  the  importance  of  doing 
the  best  possible  job  of  tracking  the  target.  In  addition,  the  simula¬ 
tion  was  both  interesting  and  challenging,  and  each  individual  wanted 
to  show  how  well  they  could  track  the  target.  It  is  therefore  assumed 
that  each  individual  was  both  well  trained  and  highly  motivated.  The 
next  step  is  to  explain  the  evaluation  process. 

4.3  Evaluation  Process 

This  section  will  outline  the  approach  for  testing  the  two  research 
hypotheses.  The  evaluation  process  involves  the  selection  of  certain 
data  parameters  and  how  these  parameters  test  the  validity  of  each 
hypothesis.  The  primary  hypothesis  is  that  the  augmentation  control 
synthesis  results  in  improved  piloted  vehicle  performance.  The  parame¬ 
ters  to  test  this  hypothesis  are  the  elevation  and  azimuth  errors  and 
the  control  surface  deflection  rates.  The  errors  will  determine  whether 
the  performance  is  better  while  the  rates  will  show  if  pilot  workload 
is  reduced.  The  parameters  compared  are  average  simulation  RMS  values 
and  analytical  RMS  values  for  both  the  unaugmented  and  augmented 
configuration.  Each  individual  flew  each  configuration  twice  and 
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the  average  RMS  values  will  be  the  average  of  all  simulations  runs 
for  that  particular  configuration.  In  addition,  pilot  comments  about 
the  task  difficulty,  task  realism,  and  any  undesirable  characteristics 
were  obtained.  These  conments  are  used  to  see  if  the  pilot  "feels"  the 
augmentation  system  is  helping  to  accomplish  the  task  better  and  to 
show  any  tendency  that  the  pilot  and  augmentation  controls  are  in 
conflict. 

The  second  hypothesis  test  involves  comparing  RMS  statistics 
between  analytical  and  simulation  results  to  determine  if  any  OCM 
parameters  need  adjustment.  The  parameters  to  be  compared  are  e^,  e^, 
Vj,  Wy,  a,  e,  XEL,  Xjg,  0,  $,  p,  q,  r,  &£,  6r,  and  Ay  .  It  is 

hoped  that  the  average  simulation  RMS  statistics  are  reasonably  close 
to  analytical  RMS  statistics.  If  they  are  close,  this  would  indicate 
the  OCM  is  indeed  doing  a  good  job  of  predicting  performance.  If  all 
the  parameters  are  grossly  off,  this  cou Id  mean  either  the  OCM  or  the 
control  synthesis  is  in  error,  or  the  simulation  is  defective.  In 
either  case,  a  closer  look  at  the  OCM,  the  control  synthesis,  and  the 
simulation  would  have  to  be  made. 

4.4  Simulation  and  RMS  Calculation  Procedures 

This  section  will  outline  the  procedures  for  operating  the  simula¬ 
tion  and  analyzing  the  data,  and  can  be  used  as  a  guide  for  other 
simulations  or  as  a  method  of  retrieving  already  stored  data.  The 
initial  step  in  the  simulation  process  is  to  insure  that  the  trim  of 
the  control  surfaces  is  set  as  close  to  zero  as  possible.  The  voltages 
of  each  control  surface  can  be  displayed  on  the  Superbee  teletype 
terminal  by  typing  $DASCAL,  and  the  voltage  values  of  various  control 
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surfaces  are  displayed.  The  first  three  channels  represent  the  eleva¬ 
tor,  aileron,  and  rudder  trim,  respectively,  and  the  appropriate  trim 
button  can  then  be  adjusted,  as  appropriate.  From  experience,  the 
voltages  should  be  less  than  .0075  volts.  The  next  step  to  actually 
start  the  simulation  sequence  by  typing  $$RWAY,  which  defines  the  data 
file,  WAYDAT,  which  contains  the  elements  of  the  system  matrix,  and 
then  starts  the  actual  simulation  program,  called  DISPLY.  A  message 
will  appear  on  the  terminal,  saying  that  the  values  of  IAUG  and  DT 
must  be  typed  on  the  terminal.  Once  these  values  are  inputted,  the 
simulation  display  wlil  appear  and  the  simulation  will  continue  until 
the  S  button  is  depressed,  which  stops  the  program  and  then  moves  the 
data,  if  any  was  taken,  to  a  working  file  location.  To  begin  copying 
data,  the  letter  C  should  be  depressed  any  time  after  the  simulation 
has  begun.  When  the  simulation  is  done,  a  message  will  appear  on  the 
terminal  screen,  which  explains  how  to  move  the  data  to  a  more  permanent 
area.  The  message  says  the  data  can  be  stored  by  typing  $C0PY,  A,  W3, 
DS(XXXXXX),  where  DS(XXXXXX)  defines  the  storage  location  and  the 
identification  scheme  is  described  in  Appendix  H. 

A  separate  program,  WAYPLT,  is  used  to  calculate  the  RMS  statistics 
for  any  parameter  and  to  plot  any  two  parameters  on  the  VERSATEC  plotter, 
if  desired.  A  more  complete  description  of  WAYPLT,  plus  a  computer 
listing,  can  be  found  in  Appendix  G.  The  first  step  is  to  select  the 
appropriate  data  file  by  typing  $SI=DS(XXXXXX).  Then,  some  output 
device  lables  must  be  set,  namely  OM  and  LO.  The  label,  OM,  should  be 
$0M=SY,  which  prints  any  output  messages  on  the  Superbee  teletype 
terminal,  while  LO,  which  determines  where  the  RMS  statistics  will  be 
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printed,  can  be  set  to  several  different  output  devices.  The  proper 
command  is  $L0SXX,  where  XX=LP  means  the  line  printer,  XX=EP  means  the 
VERSATEC  plotter,  and  XX=SY  means  the  Superbee  terminal.  Once  the 
data  file  and  output  devices  are  set,  the  program  is  called  by  typing 
$WAYPLT  on  the  terminal.  A  message  will  appear  on  the  Superbee  screen 
which  asks  for  the  values  of  IX,  IY,  IPLT,  IRMS,  and  IRATE  to  be  set. 
These  parameters  determine  which  variables  should  be  chosen,  whether 
RMS  statistics  should  be  calculated,  and  if  any  plots  should  be 
generated.  Appendix  G  gives  a  more  detailed  explanation  of  the  parame¬ 
ters,  while  Appendix  H  gives  the  variable  numbers  associated  with  IX 
and  IY.  An  example  of  the  RMS  statistics  printout  is  shown  in  Figure  5. 

4.5  Results 

This  section  contains  the  tabulation  of  results  and  a  discussion 
about  them.  Tables  14  and  15  give  the  average  RMS  statistics  of  various 
parameters  for  both  the  unaugmented  and  augmented  systems.  Table  16 
contains  a  summary  of  pilot  comments.  The  complete  tabulation  of  RMS 
statistics  and  pilot  comnents  for  each  individual  are  presented  in 
Appendix  H.  The  information  from  this  section  will  be  the  basis  for 
the  conclusions  given  in  Chapter  5. 

The  RMS  statistics  in  Table  14  are  for  three  separate,  and  differ¬ 
ent,  cases.  The  analytical  statistics  were  calculated  for  a  target  in 
a  constant  bank,  four-g  turn,  and  the  attacker  having  achieved  a 
steady-state  solution.  The  previous  simulation  case  are  the  results 
of  the  simulation  at  Wright-Patterson  Air  Force  Base,  which  was  describ¬ 
ed  in  Section  2.4,  and  while  the  data  is  based  on  a  target  in  a  constant 
four-g  turn,  the  pilot  did  not  achieve  a  steady-state  solution,  so  the 
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RMS  statistics  are  based  on  perturbations  about  some  average  value. 

The  present  simulation  results  are,  of  course,  based  on  the  target 
perturbating  about  the  constant  four-g  turn.  Because  of  these  differ¬ 
ences,  a  comparison  between  the  different  cases  should  not  result  in  a 
one-to-one  correlation. 

As  mentioned  in  Section  4.3,  the  RMS  statistics  are  an  average  of 
each  individual  flying  each  configuration  twice.  However,  for  the 
unaugmented  configuration,  only  two  of  the  four  individuals  used  the 
rudder,  and  they  had  vastly  different  rudder  deflections  and  rates. 

By  examining  the  individual  statistics,  as  given  in  Appendix  H,  the 
average  value  trends  is  the  same  as  each  individual’s  trend.  It  is 
obvious  from  the  comparing  these  averages  values  that  a  one-to-one 
correlation  does  not  exist,  but  the  differences  are  primarily  due  to 
differing  target  motions.  The  addition  of  extra  target  banking  was 
expected  to  result  in  variations  of  the  RMS  statistics,  and  the  largest 
differences  were  in  parameters  associated  with  the  attacker's  banking 
control,  such  as  e^,  idj,  <j>,  p,  r,  and  5^.  A  comparison  between  the 
unaugmented  and  augmented  configurations  shows  that  all  the  parameters 
have  a  reduction  in  average  RMS  statistics  except  Wy  and  which 
have  very  slight  increases. 

The  pilot  comments  were  based  on  the  1-10  rating  scale  mentioned 
earlier,  where  the  rating  of  10,  or  uncontrollable,  meant  the  individual 
could  not  keep  the  sight  symbol  touching  some  part  of  the  target  symbol 
throughout  the  simulation  run.  All  of  the  individuals  rated  the  simula¬ 
tion  about  the  same,  so  the  sunmary  contained  In  Table  16  provides  an 
accurate  record  of  their  comments. 
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THE  AVG  VALUE  OF  E  EL  IS-  0.03883 

THE  RMS  VALUE  OF  E  EL  IS-  0.56668 

THE  STD  DEV  OF  E  EL  IS-  0.56554 

THE  NUM  OF  POINTS-  483 

THE  AVG  'VALUE  OF  E  AZ  IS-  0.10974 

RMS  VALUE  OF  E  AZ  IS-  1.29954 

THE  STD  DEV  OF  E  AZ  IS-  1.29490 

THE  NUM  OF  POINTS-  483 

TiE  A\*G  VALUE  OF  DE  R  IS-  -0.00280 

T£  RMS  VALUE  OF  DE  R  IS-  0.45173 

'•HE  STD  DEV  OF  DE  R  IS-  0.45172 

THE  NUM  OF  POINTS-  483 

THE  AVG  VALUE  OF  DA  R  IS-  -0.00280 

THE  RMS  VALUE  OF  DA  R  IS-  0.40560 

THE  STD  DEV  OF  DA  R  IS-  0.40559 

THE  NUM  OF  POINTS-  483 

“HE  AVG  VALUE  OF  DR  R  IS-  -0.00473 

THE  RMS  VALUE  OF  DR  R  IS-  0.31750 

T-E  STD  DEV  OF  DR  R  IS-  0.31747 

THE  NUM  OF  POINTS-  483 


Figure  6 

Sample  RMS  Statistics  Printout 
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Table  14 

Unaugmented  System  Comparison 
(RMS  Performance) 


PARAMETER 

ANALYTIC 

RESULTS 

PREVIOUS 

SIM  RESULTS* 

PRESENT 

SIM  RESULTS  (AVG) 

eEL  (deg) 

1.78 

1.09 

1.33 

e/2  (deg) 

1.72 

.97 

3.80 

VT(  ft/sec) 

25.0 

34.0 

18.24 

WT  (ft/sec)! 

63.0 

48.0 

23.05 

a  (deg)  i 

1 

.97 

1.4 

.69 

6  (deg) 

.15 

.20 

.42 

*El  (deg)  j 

1.38 

1.72 

.63 

XAZ  (deg) 

.63 

2.58 

.52 

0  (deg) 

.97 

1.09 

10.03 

4>  (deg) 

1 

3.55 

4.18 

11.06 

p  (deg/sec) 

5.21 

| 

4.99 

10.02 

q  (deg/sec) 

2.58 

1.72 

1.82 

r (deg/sec) 

.57 

.92 

1.53 

«E (in) 

.23 

.27 

.20 

«A  (in) 

.15 

.22 

.52 

«R  (in) 

.10 

.22 

.20 

6e  (in/sec) 

- 

- 

.87 

6a  (in/sec) 

- 

- 

2.22 

6r  (in/sec) 

- 

- 

.88 

Ay  (ft/sec2; 

.78 

.67 

2.46 

$t(avg/rms) 

0.0/5.25 

- 

-.39/5.22 

♦Simulation  done  at  Wright-Patterson  AFB 


Table  15 


Augmented  System  Comparison 


PARAMETER 

(RMS  Performance) 

I 

ANALYTIC  RESULTS 

PRESENT 
SIM  RESULTS 

eEL(deg) 

.54 

.49 

eAZ(deg) 

.20 

1 

1.31 

Vj( ft/sec) 

4.0 

5.95 

WT( ft/sec) 

15.0 

24.58 

a  (deg) 

.36 

.57 

6  (deg) 

.06 

j 

.35 

X£L (deg) 

.33 

.72 

>AZ(deg) 

.10  i 

.17 

0  (deg) 

.17 

9.67 

<f>  (deg) 

.48 

10.03 

p( deg/sec) 

1.20 

!  2.58 

q (deg/sec) 

.97 

1.15 

r (deg/sec) 

.21 

.34 

SE(in) 

.19 

.20 

<SA(in) 

.11 

.18 

fiR  (in) 

.06 

.33 

«E (in/sec) 

- 

.36 

6A  (In/sec) 

- 

.38 

<SR  (in/sec) 

- 

.58 

Ay  (ft/sec^ 

.25 

1.87 

*t(avg/rms) 

0.0/5.25 

-.38/5.27 

Table  16 


Seminary  of  Pilot  Comments 

Was  the  augmentation  system  helpful?  All  subjects  said  Yes. 
Any  undesirable  or  annoying  characteristics?  None  of  the 
subjects  mentioned  any. 

How  would  you  rate  the  task  difficulty  from  1-10,  where 
1  =  no  inputs  required 
5  =  reasonable  difficulty  for  task 
10  =  uncontrollable 

Unaugmented?  Subjects  rated  between  5  and  7 
Augmented?  Subjects  rated  between  2  and  3 
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4.6  Sunriary 

This  chapter  has  outlined  the  evaluation  process  used  to  test  the 
two  research  hypotheses.  The  model  assumptions  were  discussed  and  it 
was  shown  that  they  were  valid  for  the  simulation.  The  actual  process 
used  to  evaluate  each  research  hypothesis  was  then  explained.  The 
procedures  for  running  the  simulation  and  calculating  the  RMS  statistics 
was  presented  next.  Finally,  the  results  of  the  simulation  were  tabu¬ 
lated  ard  discussed. 
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CHAPTER  5.  CONCLUSIONS  AND  SUGGESTIONS 


5.1  Introduction 

This  section  will  use  the  results  presented  in  Section  4.4  to  draw 
some  conclusions  about  the  research  hypotheses.  The  simulation  results 
show  that  the  control  synthesis  provides  improved  piloted  vehicle  per¬ 
formance.  Also,  comparing  the  simulation  RMS  statistics  to  analytical 
results  indicates  several  parameters  are  considerably  different.  How¬ 
ever,  the  largest  differences  are  in  parameters  associated  with  the 
attacker's  bank  control,  and  since  target  motion  was  added  to  the 
optimal  pilot  control  model  via  target  perturbation  banking,  this 
suggests  a  logical  reason  for  the  disparity.  To  test  this  premise,  the 
first  suggestion  for  future  research  is  to  do  the  simulation  with  no 
target  perturbations  and  see  if  the  RMS  statistics  yield  better  correla¬ 
tion  between  simulation  and  analytical  results.  Another  suggestion  is 
to  simulate  the  pitch  only  tracking  task  that  was  described  in  Section 
2.3. 

5.2  Conclusions 

The  primary  objective  in  this  research  is  to  determine  if  the 
augmentation  control  synthesis  results  in  improved  piloted  vehicle 
performance.  It  is  expected  that  the  RMS  statistics  will  show  a 
decrease  in  elevation  and  azimuth  errors,  since  the  methodology  requires 
it.  However,  it  is  hoped  that  pilot  workload  will  also  be  decreased. 


and  the  pilot  workload  was  measured  objectively  by  using  RMS  statistics 
on  the  control  surface  deflection  rates*  as  well  as  subjectively,  based 
on  pilot  cormients.  The  average  RMS  statistics  for  both  the  augmented 
and  unaugmented  configuration  are  shown  in  Table  16,  and  the  elevation 
and  azimuth  errors  are  both  reduced,  especially  in  azimuth.  The 
deflection  rates  for  the  aileron  and  elevator  are  reduced,  as  well  as 
for  the  rudder.  The  lower  deflection  rates  show  a  reduced  piloted 
workload  and  this  premise  is  supported  by  the  favorable  pilot  comments, 
which  were  summarized  in  Table  15.  All  the  pilots  rated  the  task 
difficulty  significantly  easier  in  the  augmented  case  and  they  all  felt 
that  there  were  no  unfavorable  characteristics  for  the  augmented  system. 
Because  of  the  lower  RMS  statistics  and  favorable  pilot  corments,  it  is 
felt  that  the  piloted  vehicle  performance  is  improved  significantly  by 
the  augmentation  system. 

The  secondary  research  objective  was  to  compare  the  analytical  RMS 
statistics  of  the  OCM  to  simulation  statistics  in  order  to  determine  if 
any  significant  parameter  changes  are  required.  The  analytic  results 
are  based  on  a  constant  target  bank  angle,  so  it  is  not  surprising  that 
the  attacker  bank  parameters  do  not  match  very  well.  Also,  note  that 
Xpj,,  the  sight  azimuth  angle,  shows  excellent  agreement,  and  since 

eAz =  xaz  _baz’  the  d1sPanty  in  eaz  must  ^  Priman1y  due  to  b^»  the 

target  azimuth  angle,  and  this  indicates  the  target  bank  is  largely 
responsible  for  the  disparity.  It  is  therefore  felt  that  any  deci¬ 
sion  about  OCM  parameter  adjustment  should  use  a  similar  analytical  study. 
Some  observations,  however,  can  still  be  made  about  the  parameters. 

It  should  be  noted  that  the  parameters  associated  with  aileron  and 
rudder  deflections  and  deflection  rates  are  considerably  higher  than 


the  elevator  parameters,  which  indicates  the  pilots  have  a  much  more 
difficult  time  of  controlling  the  azimuth  error.  Also  the  sideslip 
angle,  6,  and  lateral  acceleration.  Ay.  are  clues  that  would  be  avail¬ 
able  in  a  motion  simulator,  and  since  this  simulation  is  fixed  based, 
it  is  not  surprising  that  the  simulation  results  are  much  larger  than 
analytic  results.  Finally,  pilot  comments  from  DJB  and  WAY,  both  who 
have  air-to-air  tracking  experience,  indicate  the  simulation  is  fairly 
realistic  of  actual  multi  axis  tracking  tasks.  It  is  therefore  felt 
that  the  simulation  is  realistic,  and  while  the  RMS  statistics  do  not 
match  exactly,  the  introduction  of  the  target  perturbation  angle 
accounts  for  the  disparity. 

5.3  Suggestions 

This  section  contains  three  suggestions  for  possible  follow-on 
research.  The  first  suggestion  is  fairly  obvious  and  relatively  simple 
to  implement.  An  analytical  study  with  the  additional  target  motion 
added  could  be  done  and  the  calculated  RMS  statistics  could  then  be 
compared  to  the  present  simulation  results.  A  second  suggestion  is  to 
simulate  the  pitch-only  tracking  task  described  in  Section  2.3.  This 
type  of  simulation  could  be  used  to  determine  if  a  pilot’s  strategy  for 
correcting  elevation  would  be  the  same  for  both  a  pitch-only  and  a 
multiaxis  tracking  task.  A  third  possibility  would  be  to  use  some 
more  advanced  simulation  facility,  such  as  the  Flight  Dynamics  Labora¬ 


tory  at  Wright-Patterson  AFB  to  incorporate  motion  cues. 
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Table  17 

Unaugmented  vs.  Augmented  System  Comparison 
(RMS  Performance) 


PARAMETER 

UNAUGMENTED 

AUGMENTED 

eEL 

1.33 

.49 

eAZ 

3.80 

1.31 

*E 

.87 

.36 

SA 

2.22 

.38 

.88 

.58 

I 

if 


5.4  Summary 

This  chapter  seems  to  validate  the  decision  to  use  the  optimal 
pilot  control  model  for  modeling  the  multi  axis  air-to-air  tracking  task. 
The  simulation  results,  as  shown  by  a  reduction  in  elevation  and  azimuth 
errors,  as  well  as  reduced  control  surface  deflection  rates,  indicate 
the  piloted  vehicle  performance  characteristics  are  significantly  improv¬ 
ed  by  using  the  augmentation  control  synthesis.  Pilot  comments  about 
the  piloted  vehicle  characteristics  also  support  this  idea.  While  the 
RMS  statistics  do  not  match  the  analytic  statistics,  the  introduction 
of  target  perturbation  motion  accounts  for  the  disparities.  One 
suggestion  for  future  research  is  to  delete  the  target  motion  and  then 
compare  simulation  results  to  analytic  results.  Other  suggestions 
include  simulating  the  pitch  tracking  task  and  simulating  the  multiaxis 
task  on  a  more  sophisticat'd  simulation  system. 

It  should  be  remembered  that  this  research  is  not  an  end  in  itself, 
but  only  one  small  step  in  a  larger  plan.  The  ultimate  goal  is  to  find 
some  theoretical  method  of  accurately  predicting  piloted  vehicle  perfor¬ 
mance,  especially  when  higher  order  system  dynamics  are  present.  Since 
classical  methods  are  inadequate  for  accurately  predicting  the  piloted 
vehicle  performance  in  this  situation,  the  optimal  control  model  was 
chosen  as  one  method  for  predicting  the  performance  characteristics. 

An  augmentation  control  synthesis  method  was  then  employed  to  achieve 
improved  analytic  performance,  but  additional  simulation  results  were 
required  to  validate  this  hypothesis.  The  primary  goal  of  this  research 
was  to  simulate  a  multi  axis  air-to-air  tracking  task,  whose  analytical 
development  had  previously  been  done  by  Schmidt. A  secondary 
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objective  was  to  use  these  simulation  results  to  determine  If  any 
changes  to  the  OCM  parameters  were  required  to  better  model  the  task. 

The  simulation  RMS  statistics,  as  well  as  pilot  comments,  indicate  the 
augmentation  synthesis  does  indeed  provide  improved  piloted  vehicle 
performance.  In  addition,  although  the  simulation  statistics  do  not 
closely  match  all  the  analytic  statistics,  the  addition  of  target 
perturbation  motion  explains  the  differences.  Therefore,  the  optimal 
pilot  control  model  seems  one  promising  method  for  accurately  predicting 
piloted  vehicle  performance  characteristics  and  the  augmentation  control 
synthesis  provides  one  method  for  improving  the  piloted  vehicle  perfor¬ 
mance  . 


i  1 1  nJ  iill'IiH 
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APPENDIX  A 

State  System  Variable  Matrices 


This  appendix  contains  printouts  of  all  the  appropriate  matrices. 
The  A,  B,  C,  and  D  matrices  are  the  linear  system  matrices  defined  by 


x  =  Ax  +  Bu 


y(t+At)  =  Cx  +  Du 


The  Kx  and  Ku  matrices  are  the  gain  matrices  calculated  by  the  control 
synthesis,  and  the  AAUG  and  BAUG  matrices  were  calculated  using 


AAUG  =  A  -  BK 


BAUG  =  B  -  BKu  . 


The  EXP  TAU  matrix  is 


‘'At, 
e  dT 


and  EXP  TAU  times  B  matrix  is 


(fVdx)B 

*  O 


The  transition  matrix  is  e  . 
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THE  KX  AUGMENTATION  MATRIX 


-a. 0330000 

-0.0050000 

-0.0030000 

-0.0004000 

0.0090000 

0.0050000 

-C. 0048000 

-0.0000400 

-0.0000200 

-5.3109998 

0 . 1270000 

0.0750000 

1480000 

5.4699997 

3.5849997 

8.2230300 

-0. 1050000 

-0.0810000 

-S . 1 170900 

-1.9770000 

-0.3340000 

5.3119997 

-0.0310000 

-0.0260000 

1370C00 

0.6440000 

-0.0010000 

t .3430000 

-9.4649996 

7.5169997 

0.0230000 

-1.8380001 

-0.0479000 

0.0330000 

-2.2059998 

-5.7009992 

THE 

KU  AUGMENTATION  MATRIX 

0.3180000 

-0.0020000 

0. 
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APPENDIX  B 

State  Transition  Matrices  Calculations 


.AAt 


Two  state  transition  matrices,  e  and  (j 
puted  to  update  the  x- vector.  The  series  representation  of  enT  is  given 
by 


eAx)B,  must  be  precom- 
At 


At  —  T  4.  A  4.1  *22.1  .3  3. 

e  =  I  +  At  +  ^  A  t  +  ti  A  T  + . 


I! 


(Al-1) 


Integrating  term  by  term,  over  the  limits  0  to  At,  yields 

|  eArdt  =  iAt  +  -^-j  A(At)2  +  -j|  A2(&t)3  +  ...  .  (Al-2) 

Note  that  multiplying  (Al-2)  by  the  matrix.  A,  and  adding  the  identity 
matrix,  I,  gives 


I  +  AAt  +  |,  A2(At}2  +  1  A3(At)3  +  . 


(Al-3) 


.AAt 


that  must  actually  be  computed  is 


AAt 


which  is  the  series  representation  of  e  .  Therefore,  the  only  matrix 

fAt 

dr.  Then  multiply  it  by  A  and 
°r&t  At 

e  dt  by  the  B  matrix  to  obtain 
'  o 

the  second  state  transition  matrix. 

e  dt,  was  computed  by  adding  additional  terms  of 


add  I  to  get  e  ,  and  multiply 


The  matrix, 

the  series  until  the  contribution  from  the  higher  order  terms  was  negli¬ 
gible.  The  first  term  in  the  series  was  called  El,  and  the  second  term, 
E2.  If  the  difference  of  every  element  of  the  resultant  matrix  is  less 


r 

in  m*""  ■  -  — -  r  in  hr  it  T  -iny-T 


a  specified  value,  .000001  in  this  program,  then  the  lower  term  is 
r At  . 

added  to  eATdx,  and  the  computation  is  complete.  If  any  element  is 

0  (At  a 

greater  than  the  specified  amount,  the  lower  term  is  added  to  I  eMTdx 

••o 

and  the  matrix  El  is  set  equal  to  E2.  The  next  higher  term  is  then 

computed  by  multiplying  the  new  El  matrix  by  AAt  and  then  dividing  by 

the  order  of  the  term.  For  example,  the  third  term  in  the  series  is 

obtained  by  multiplying  the  second  term  by  ~  AAt.  This  higher  term  is 

now  the  new  E2  matrix  and  the  procedure  outlined  above  is  repeated  until 

the  higher  term  produces  no  appreciable  change.  The  accuracy  of  these 

calculations  is  limited  only  by  computer  accuracy,  and  since  they  are 

precomputed  matrices,  they  can  be  calculated  on  any  computer,  to 

achieve  the  desired  level  of  accuracy.  As  mentioned  earlier,  the  two 

/■At  . 

required  state  transitions  are  now  obtained  by  multiplying  eAtd t 

/■At 

by  A  and  adding  the  identity  matrix  and  by  multiplying  eHTdt  by  the 

•*0 

B  matrix. 
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APPENDIX  C 

Sum  of  Sinusoids  Control  Inputs 


The  equation  used  to  obtain  the  target's  perturbation  aileron 
deflection  is  given  by 

13 

6A(t)  =  l  A.  sin  (w.t+A.)  ( Cl - 1 ) 

i=l  1  11 


where  the  values  of  A^ ,  w. ,  arid  <j>.  are  shown  below. 


2 

1 

0.001 

2 

0.001 

3 

0.001 

4 

0.001 

5 

0.05 

6 

0.05 

7 

0.1 

8 

0.1 

9 

0.1 

10 

0.1 

11 

0.1 

12 

0.1 

o)i  (rad/sec) 

.18850 
.31416 
.50265 
.87965 
1.44513 
2.31628 
3.07876 
4.20973 
5.78053 
8.23097 
11.2469 
15.77079 
23. 93894 


<t>.j  (rad) 
ir/6.5 

2<|>i 

3<t>.j 

4<f>-| 

5(j>i 

6<t>^ 

7*1 

8*i 

9<|>i 

io*i 

n*-| 

12*-| 

13*1 


#  of  cycles  in  100  sec.  run 

3 

5 

8 

14 

23 

34 

49 

67 

92 

131 

179 

251 


13  0.1 


381 
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The  values  for  u^,  the  frequencies,  were  selected  so  that  no  frequency 
was  an  integer  multiple  of  any  of  the  other  twelve  frequencies.  The 
magnitudes.  A..,  were  selected  by  experimentation  to  give  an  rms  value 
of  4>y,  using  the  target  motion  equations,  of  approximately  5  degrees. 
The  magnitudes  selected  above  resulted  in  an  rms  value  of  <j>T  =  5.25 
for  a  At  =  .0887  and  a  total  time  of  90  seconds. 
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APPENDIX  D 

Coordinate  Transfo rmation  Calculations 


The  transformation  from  the  target's  coordinate  system  (TA,  JA»  RA), 
to  the  attacker's  coordinate  system  (?A,  R^),  is  obtained  using  the 

transformation  matrix  relationships  of 

0,  87) 

(Bl-1) 


where  the  individual  transformation  matrices  are 


f  COS  0 

0 

sin  0 

_COS  t|j 

-sin 

0~ 

T(0)  =  |  0 

1 

0 

TU)  = 

sin  ii> 

cos 

* 

0 

j  -sin  0 

0 

COS  0 

0 

0 

1_ 

!  i  o  o 

i 

T(<|.)  =  I  0  cos  <t>  -sin  $ 

» 

i 

|  0  sin  <(>  cos  $ 

The  resultant  transformation  matrix  is  represented  by 


(Bl-2) 


(Bl-3) 
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One  simplification  is  obtained  by  remembering  that  T(-tf,A)  T(^T)  = 

T(-(i|»A-i|>y))  =  T(-Ai^) .  Also,  the  value  of  ^  is  equal  to  4>^  +  4>|, 
where  4>1  =  the  steady  state  target  bank  angle,  which  is  75.5°,  while 
<(>j  is  the  target's  perturbation  bank  angle.  Substituting  the  values 
into  the  appropriate  matrices,  and  using  a  small  angle  approximation 
for  yields: 

Til  =  cos(eA)  cos ( Aip)  cos(Qj.)  +  sin(0A)  sin  (ej 

T12  =  cos(©A)  cosU^sIntej-Hsinf^)  +  cos^)]  +  cos(eA)  sin^). 
[cosf^)  -  4>|  sin($-j)]  -  sin(0A  )  coi(&j.)[sin(4^ )  +  ^  cos(  $.)] 

T13  =  cos(0A)  cos ( Aip)  sinte^J^cos ( ^ )  -  *j.  sin^)]  -  cos(0A)  . 

sin( a^) [sin( ^ )  -  <^cos($1)]  -  sin(0A)  cos( aj.)[cos( ^ )  -^sinf^  )] 

Tgf  =  cos(Oj.)[sin($A)  sin(0A)  cos(ai^)  -  cos ( <t»A)  sin ( A\p) ]  -  sin(<|>A)cos(oA) 
sinfey) 

T22  =  sin(0y)[sin(<|>A)  sin(0A)  cos(A<|>)  -  cos(<j>A)  si n ( A\p) ] [sin ( )  + 
(j'yCost^)]  +  [ s i n ( 4>a )  sin(©A)  sin(A^)  +  cos(<f>A)  cos(a^)]  . 

[cosf^)  -  s i n ( 4>1 ) ]  +  sin($A)  cos(0A)  cos(oT)]sin(4>1 )  + 

<t>j  COS4>^] 

T23  =  sin(©T)[sin( <j>A)  sin(©A)  cos(a^)  -  cosUA)  sin(A*)][cos(4>1 )  - 
4>y  sin^)]  -  [si n ( «t>A)  sin(0A)  sin(A^)  +  cos(<j>A)  cos(a^)]  . 

[sin(^)  +  <ji|  cos  (*1  )3  +  sinf^)  cos(0A)  cos  (Gy)  [cos  )  - 
sin($1 )] 
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T31  =  cos(©T)[cos(<|>A)  sin(e^)  cos(a^)  +  si n ( 4>A)  sin(A^)]  - 
cos(4>a)  cos(qa)  sin(eT) 

T32  =  sin(©j)[cos(<^A)  sin(©A)  cos (a^)  +  si n (4>A)  sin(Aij<)][sin(4>1 )  + 
fj.  cos^)]  +  [sin ( 4>a)  sin(©A)  sin(A^)  +  cosU >A)  cos(a^)]  . 

[cos ( 4>i )  -  fj-  sin ( 4>1 )]  +  cos ( 4>A)  cos(eA)  cosfeyKsin^ )  + 
cosf^)] 

T33  =  sin(eT)[cos(4iA)  sin(eA)  cos(a^)  +  si n ( <t>A)  sin(Ag>)][cos($.j )  - 
sin($1)]  -  [cos(4>a)  sin(©A)  sin ( A^)  -  sin(<t>A)  cos(a^)]  . 

[sin($1)  +  <)>j  cos^)]  +  cos ( 4>A)  cos(©A)  cos(©t)[cos(<!>1  )  -  c^si n ( ) ] 

The  primary  use  of  the  transformation  is  to  transform  vT  and  wT  , 

't  T 

the  perturbation  velocities  of  the  target,  in  the  target's  coordinate 
system,  into  vT  1  and  wT  ',  the  same  quantities,  but  in  the  attacker's 

iT  iT 

coordinates.  The  equations  for  vT  '  and  wT  '  are  given  by 

*T  *T 

Vj  '  =  [g  cos ( 4>i ) ]  <t>T’  jt  =  8.05  <(>T'  (Bl-4) 

wr  '  =  [-g  sin  ( ) 3  <t>j'  k  =  -31.18  <j>T'  (Bl-5) 

The  transformation  is  found  by  using 


APPENDIX  E 


Graphics  Subroutines  and  Entities 

The  FORTRAN  subroutines  used  to  generate  the  simulation  display, 
with  a  brief  explanation  of  each,  are  presented  below.  A  complete 
description,  plus  many  other  grahpics  instructions,  can  be  found  in 
Reference  4. 

1 .  GSPIN 

GSPIN(ICA,  NWDS,  IER)  establishes  an  array,  called  ICA,  in  a 
permanent  core  memory  location.  The  dimension  of  ICA  is  set  by  NWDS, 
and  ICA  will  provide  the  area  for  creating,  displaying  and  manipulating 
the  various  entities.  IER  is  an  error  return  variable. 

2.  GDRAW 

GDRAW(LX,LY,  NPTS,  LCODE)  creates  the  code  required  to  draw  one 
or  more  straight  lines  between  a  set  of  points  contained  in  the  arrays 
LX  and  LY.  NPTS  determines  the  number  of  points  and  LCODE  is  a  four 
digit  integer  that  determines  the  exact  display.  For  example,  LCODE  = 
4000  tells  the  vector  generator  to  draw  line  segments  between  the 
points,  and  the  values  in  the  Lx  and  LY  arrays  correspond  to  raster 
units  on  the  display  screen. 

3.  GVECT 

GVECT(IVECT,  X,  Y)  creates  and  updates  the  vector  entity  numbered 
IVECT.  Lines  are  drawn  from  the  point  (X(l),  Y( 1 ) )  to  (X(2),  Y(2) )  to 


I 


* 
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(X(3),Y(3))  and  so  on  for  all  the  points  contained  in  the  storage 
arrays,  X  and  Y. 

4.  GSTOP 


As  might  be  expected,  this  subroutine  shuts  off  the  simulation 
display  to  prevent  screen  burnout. 


The  graphics  entities  used  in  the  programing  are 


ENTITY  NUMBER 
1 
2 

3 

4 

5 


CONTENTS 
Sight  Symbol 
Target  Symbol 
Upper  Horizon  Line 
Middle  Horizon  Line 
Lower  Horizon  Line 
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APPENDIX  F 

Simulation  Computer  Program 


This  appendix  contains  the  FORTRAN  computer  program,  called  DISPLY, 
that  simulates  the  multiaxis  tracking  task.  First,  the  computer  vari¬ 
ables  will  be  listed  and  defined,  then  the  subroutines  called  in  the 
program  will  be  given  with  a  brief  description  of  each,  and  finally,  a 
computer  listing  of  DISPLY  will  be  given. 


F.1.1  Computer  Variables 
SYMBOL 


MEANING 


A,  AAUG 
AA 

ADT 

AGN 

AK 

B,  BAUG 


Unaugmented  and  augmented  system  matrices,  A  and  A=11„ 

13  au9 

Array  of  A^ ,  used  to  calculate  6.  =  T  A.sin(w.t+<j>. ) 

i  «  i=l  1  1  1 

The  matrix,  A,  multiplied  by  at 


aug 


Aileron  gain  factor 

Weighting  factor  in  Runga-Kutta  integration  routine 
Unaugmented  and  augmented  system  matrices,  B  and  B 
BKx,  BKu,  Bx,  Bu  The  B  matrix  multiplied  by  the  Kx  and  Ku  matrices,  and 
the  x  and  u  vectors,  respectively. 

C  The  system  matrix,  C. 

CK  Weighting  factor  in  Runga-Kutta  integration  routine 

CPYFLG  Logical  variable  used  to  decide  if  data  will  be  taken 

CX  The  C  matrix  multi  pled  by  the  X- vector 

D  The  system  matrix,  D 

DAL  The  target  aileron  deflection,  in  radians 
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DT 

DTT 

DX,  DU 

El,  E2 

EAT,  EATX 

ERRT 

ETAU 

ETAUB,  ETAUBU 
HX.HXl ,HX2,HX3 
HY  ,HY1  ,HY2  ,HY3 
II,  12 

I  AUG 
IMAGE 
KU,  KX 
KUU 

KXU,  KXX 
OMEGA 

PHIREF.PHIREL 

POSX, P0SX2,P0SX3 

POSY. POSY2 ,POSY3 
RAD  755 

RDCF 


The  iteration  step  size.  At 

Array  of  the  step  sizes  used  in  the  Runga-Kutta 

integration  routine 

The  D  matrix  multiplied  by  the  X  and  Y  vectors, 
respectively 

Storage  matrices  used  to  calculate  state  transition 
matrices 

e^At  and  e^^  times  the  X  vector,  respectively 
Array  of  errors,  used  in  Runga-Kutta  routine 

iv* 

j  0 

ETAU  multiplied  by  the  B  and  BU  matrices,  respectively 
Storage  arrays  for  horizon  line  endpoints  (x-di recti  on) 
Storage  arrays  for  horizon  line  endpoints  (y-di recti  on) 
Logical  variables  used  to  determine  if  printing  should 
be  done 

Logical  variable  used  to  decide  if  augmentation  used 

Storage  array  for  graphics  instructions 

The  augmentation  matrices,  K  and  K  ,  respectively 

The  KU  matrix  times  the  U- vector 

The  KX  matrix  times  U  and  X-vectors,  respectively 

13 

Array  of  w. ,  used  in  6.  =  l  A.  sin(w.t+<|). ) 

i  *  i«l  1  11 

,  the  attacker’s  and  target  bank  angles 

Arrays  for  storage  of  vector  endpoints  (x-direction) 

Arrays  for  storage  of  vector  endpoints  (y-direction) 

75.5°,  converted  to  radians 

Conversion  factor,  radians  to  degrees 


AD- A 125  302  COMPUTER  SIMULATION  OF  A  MULTIAXIS  AIR-TO-AIR  TRACKING  Id 
TASK  USING  THE  OPT.\  (U)  AIR  FORCE  INST  OF  TECH  ^ 

WRIGHT-PATTERSON  AFB  OH  W  A  YUCUIS  DEC  82 
UNCLASSIFIED  AF I T/CI /NR-82-62T 


F/G  9/2 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 
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REF.REFSX.REFSY 

RVALS 

RX.RY 

SX.SY 

TAUR 

TGT 

TGTX.TGTY 

TX.TY 

U 

UAUG 

X 

XBUF 

Y 


Arrays  for  reference  positions  on  the  display  screen, 
in  raster  units 

Array  of  pilot  control  inputs,  in  volts 
Array  of  vector  endpoints  of  reference  cross 
Array  of  vector  endpoints  of  sight  symbol 
Target  aircraft  roll  constant,  tr 
Array  of  target  motion  variables,  6^,  pT,  ^ 
Array  of  vector  endpoints  of  target  symbol,  with 
bank  included 

Array  of  vector  endpoints  of  target  symbol,  with 
no  banking 

Array  of  attacking  pilot's  control  inputs,  up 
Array  of  equivalent  augmented  control  inputs,  ua 

fl 

Array  of  state  vector,  x 

Storage  matrix  for  data 

Array  of  the  observation  vector,  y 


F.2.  Simulation  Subroutines 

a.  DAS IN  -  converts  analog  voltage  readings  of  strain  gages, 
used  to  calculate  stick  forces,  into  a  digital  voltage,  by 
averaging  four  readings  over  a  1  msec  time  period 

b.  FREMAT  -  allows  for  free  format  Inputting  of  variables 

c.  FSDIO  -  transfers  a  data  storage  matrix,  XBUF,  to  a  buffer 
location,  and  then  allows  the  disk  to  have  direct  access 
to  the  data 

d.  GMPRD,  GMADD,  GMSUB  -  multiplies,  adds,  and  subtracts  two 
matrices,  respectively 
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e.  TSECS  -  reads  out  the  real  time 

f.  TSET  -  starts  the  real  time  clock 

F.3.  Computer  Listing 


F 


C**  DISPLY  -  GIVES  SIGHT /TARGET  DISPLAY 
DIMENSION  IMAGE( 1500) fRVALS(4> 

LOGICAL  A0RDC.CPYFLG 

REAL  flK< 12) ,BU( 12) ,X( 12) ,Y( 10) , AAUGt 12, 12) 

REAL  BAUGC12,3),U<3),XBUF(28,7) 

REAL  EAT<12,12) ,ETAU( 12,12) |ETAUB(12|3) , El C 12,12) ,E2(12,12) 
REAL  ADT( 1 2 , 1 2) , EATXC 1 2) , ETAUBUC 1 2) 

REAL  P0SX3C  6) , P0SY3( 6) 

REAL  A(12,12),IO<X(3),KUU(3),UAUG(3) 

REAL  B(12(3) ,KX(3,12) ,KU(3,3) 

REAL  BKX( 12, 12) ,BKU( 12,3) , C( 10, 12) ,CX( 10) ,D( 10,3) 

REAL  DU( 10) ,TGT(5) ,DTT(5) ,ERRT(5) 

REAL  RX(4) ,RY(4) ,TX(S) ,TY(S) ,REF(4) 

REAL  SX( 8)  ,SY(  8) , POSXC  B) , POSY(B) , P0SX2C 2) , P0SY2( 2) 

REAL  HX(2) ,HY(2) ,HX1<2) ,HX2(2) ,HX3(2) ,HY1 (2) ,HY2(2) ,HY3(2) 

REAL  AK(4),BK(4),CK(4) 

REAL  AA(13) ,OMEGA( 13) ,PHIT2( 13) 

I NTEGER  SI , DO , ST ACK ( 1 0 ) , LO , SH I  FT , C ASE 
DATA  AA/4*0 .001,3*0.05,6*0.1/ 

DATA  OMEGA/ . 1 885 , . 3141 6 , . 50265 , . 87965 ,1.44513,2.31628,3. 07876 , 
C  4 . 20973 , 5 . 78053 , 8 . 23097 , 1 1 . 2469 , 1 5 . 77079 , 23 . 93894/ 

DATA  NCHNLS/X ' F000 ' / 

DATA  SI/1/, LO/5/, SHIFT/8192/, DO/50/ 

DATA  STACK/14, 1,1, 7*0/ 

DATA  RX/460. 0,560. 0,510. 0,510.0/ 

DATA  RY/765. 0,765. 0,715. 0,815.0/ 

DATA  TX/-80. 0,80. 0,-20. 0,0. 0,0. 0,20.0/ 

DATA  TY/0. 0,0. 0,0. 0,30. 0,30. 0,0.0/ 

DATA  SX/-20. ,-5. ,0. ,5. ,20. ,20. ,-20. ,-20./ 

DATA  SY/20. ,20. ,30. ,20. ,20. ,-20. ,-20. ,20./ 

DATA  HX/-400. 0,400.0/ 

DATA  HY/0. 0,0.0/ 

DATA  REF/B10. 0,510. 0,210. 0,765.0/ 

DATA  AK/.5, .292892319, 1 .70710678, .166666667/ 

DATA  BK/2.,1.,1.,2./ 

CKU)-AK(l) 

CK(2)-AK(2) 

CK(3)-AK(3) 

CK(4)« .5 

REFSX-REFC 2)-300 .0*0.031 

REFSY-REF(4)-300 . 0* . 1 27 

REFTX-REFSX 

REFTY-REFSY 

PI -3. 141 592654 

N-0 

RNUM-0.0 

250  RNUM-RHUmi  .0 
N-N+l 

PHIT2(  N) -RNUmPI /6 . 5 
IFCN.LT.13)  GO  TO  250 

C  DEFINE  DATA  FILE  AND  SET  UP  IMPUT  DEVICE 

DEFINE  FILE  DO(5000,12,U,IREC) 

CALL  A*I NT ( STACK ,'SY') 

ISY-IPUW1('SY')+32 

CASE-'A'+SHIFT 


i 

i 
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C  SET  CONCURRENT  DISK  WO  PARAMETERS 

ISECT-0 
INDEX-0 
IURITE-2 
IMMED-1 
I READ- 1 
IUAIT-0 

C  ZEROING  OUT  THE  MATRICES 

DO  10  1-1,12 
DO  11  J-1,12 
El(I,J)-0.0 
ETAU( I ,J)-0.0 
A( I , J)-0.0 

IF( J.LT.4)  KXtJ,I)«0.0 
IFt J.LT .4)  U( J)-0.0 
IFt J.LT.4)  BKU( 1 , J)*0.0 
BKX( I , J)-0.0 
IFCI.LT.il)  CXtI)-0.0 
X(I)-0.0 

IF(I.LT.ll)  Yt I )-0.0 
IF< J.LT.4)  B< 1 , J)-0.0 
IF(J.LT.ll)  C( J,l)-0.0 

11  CONTINUE 
10  CONTINUE 

DO  12  1-1,10 
DO  13  J-1,3 
IF< I .LT .4)  KUt I , J)-0.0 
IFtI.LT.6)  TGT( I )-0.0 
I F ( I . LT . 6)  DTTt  I )-0.0 
IF<I .LT.6)  ERRTt 1 )«0,0 
13  D( I , J)-0.0 

12  CONTINUE 

C  IMPUT 

READtSI , 1 ) 

READ(SI.l) 

READ(SI.l) 

READ(SI.l) 

READ(SI.l) 

READ! SI ,1 ) 

READtSI , 1 ) 

READ(SI.l) 

READtSI ,1) 

READ(SI.l) 

READtSI ,1) 

READtSI ,1) 

READt  SI , 1 ) 

READtSI ,1) 

READtSI ,1) 

READt SI, 1) 

READtSI, 1) 

READt SI, 1) 

READt SI, 1) 

READtSI ,1) 

READtSI,!) 


LINEARIZED  AIRCRAFT^SIGHT  DYNAMICS 

At  1 ,9) ,At 1 , 12) ,At2,6) ,At2,7) ,At2, 1 1 ) ,At2, 12) 

At  3,6) ,At3,7) , At3,9) ,At3,ll) ,At4,3) ,At4,4) 

At4,5) ,At4,8) , At4,9) ,At4, 1 1 ) , At 4, 12) 

At5,2) , At5,3) , At 5,4) ,At5,5) ,At5,B) , Af5,10) 

At5,l l) ,At5,12) ,At6,6) ,At6,7) ,At6,9) ,At6,12) 

At  7, 6), At  7, 11), At 8,7), At  8, 8), AtB, 9) 

At9,7) ,At9,8) ,At9,9) ,At9,l 1 ) ,At9,12) 

At  10, 7), At  10, 10), At  10, 11), At  10, 12), At  11, 9), At  11, 10) 
At  11, 11), At  11, 12), At  12, 9), At  12, 10), At  12, 11), At  12, 12) 
Bt8,l ) ,Bt 9,1 ) ,Bt 10,2) ,Bt 10,3) 

Bt 1 1 ,2) ,Bt 1 1 ,3) ,Bt 12,2) ,Bt 12,3) 

Ct 1 , 1 ) ,Ct 1 ,3) ,Ct 1 ,4) ,Ct 1 ,6) ,Ct 1 ,6) 

Ct2,3) ,Ct2,4) ,Ct 2,5) ,C(2,6) ,Ct 2,7) ,Ct2,8) 

Ct 2,9) ,Ct2, 1 1 ) ,Ct 2, 12) ,Ct3, 1 ) ,Ct3,2) ,Ct3,3) 

Ct 3,5) ,Ct3,6) ,Ct 3,7) ,Ct3,8) ,Ct 3, 10) 

Ct4,2> ,Ct4,3) ,Ct4,4) ,Cf4,5) ,Cf4,6) ,Ct4,7) 

Ct4,B) ,Ct4,9) ,Ct4, 10) ,Ct4, 1 1 ) ,Ct4,12) 

Ct  5 , 1 ) ,Ct5,3) ,CtS,4) ,Ct5,6) ,Ct5,8) 

Ct6,6) ,Ct6,7) ,Ct6,B) ,Ct6,9) ,Ct 6, 1 1 ) 

Ct7, 1 ) ,Ct7,2) ,Ct7,3) ,Ct7,5) 

Ct7,6) ,Ct 7,7) ,Ct 7, B) ,Ct 7, 10) ,Ct 8,6) ,Cf 8,7) 


'  e' 


•  “•  "S'  .*  .. 
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READCSI, 1)  CC8,8),CC8,9),CC8,10),CCB,11),CC8,12) 

READCSI.1)  CC9,9) ,CC9,10) ,CC9, 1 1 ) ,CC9, 12) ,CC 10,7) 

READCSI , 1 )  DC2,1),DC4,1),DC4,2),DC4,3) 

RE ADC 51 , 1 )  DC 6,1), DC 8,1), DC 8, 2), DC 8, 3) 

READC5I , 1 )  KXC 1,1) ,KXC 1,2) ,KXC 1 ,3) ,KXC 1 ,4) ,KXC 1 ,5) ,KXC 1 ,6) 
READC SI , 1 )  KXC 1 ,7) ,KXC 1 ,8) , KXC 1 ,9) ,KXC 1 , 1 B) , KXC 1 , 1 1 ) , KXC 1 , 12) 
READCSI ,1)  KXC2.1) ,KXC2,2) ,KXC2,3) ,KXC2,4) ,KXC2,5) ,KXC2,6) 
READC SI , 1 )  KXC 2, 7) , KXC 2, 8) , KXC 2,9) , KXC 2,18) , KXC 2,11) , KXC 2, 12) 
READCSI, 1)  I0<C3,1 ) , KXC 3, 2) , KXC 3,3) , KXC 3, 4) , KXC 3, 5) , KXC 3, 6) 
READCSI , 1 )  KXC 3, 7) , KXC 3, 8) , KXC 3 ,9) , KXC 3, 18) , KXC 3, 1 1 ) ,KXC3, 12) 
READCSI , 1 )  KUC 1,1) ,KUC 1 ,2) ,KUC2, 1 ) ,KUC2,2) ,KUC3, 1 ) ,KliC3,2) 
READC SI, 1)  ESN , AGN , RGN , AUG 
READC SI, 2)  II 

1  F0RMATC8F18.4) 

2  F0PHATC8I 18) 

C  SET  UP  AND  INITIALIZE  DISPLAY  AREA 

CALL  GSPINC IMAGE, 1588, I ER) 

CALL  GDRAUCRX,RY,4,4800) 

3  CALL  FREMAT 

C  READ  IN  IF  AUGMENTATION  IS  REQUIRED 

C  CSET  I AUG-0.0  IF  NO  AUGMENTATION  REQUIRED) 


C  READ  IN  DT,  THE  STEP  INTERVAL  SIZE 

URITEC ISY.4) 

4  FORMAT  C  '*0°********^^ 

C  '***  ENTER  AUG  OPTION, DT  ***',✓, 

C  '***  AUG  OPTION-8  MEANS  NO  AUG  now.',/, 

C  '***  AUG  OPTION-1  MEANS  LVL  2  AUG  ***',/, 

C  '***  DT-STEP  SIZE-SHOULD  BE  >  0.8807  ***',✓, 


READC1SY,5)  IAUG.DT 
5  FORMATCV) 

RNUM2- 1 . 0 

C  IF  I AUG-0  NO  AUGMENTATION  REQUIRED  AND 

C  AAUG.BAUG  ARE  SET  EQUAL  TO  A,B  MATRICIES 

IFC IAUG.NE.0)  GO  TO  14 
DO  14  1-1,12 
DO  15  J-1,12 
AAUGC I , J)-AC I , J) 

IFC J.LT.4)  BAUGC I , J)*BC I , J) 

15  CONTINUE 
14  CONTINUE 

C  IF  I AUG  N.E.  0,  AUGMENTATION  IS  REQUIRED  AND 

C  AAUG  AND  BAUG  ARE  CALCULATED 

IFC IAUG.EQ.0)  GO  TO  16 
CALL  GMPRDCB,KX,BKX,12,3,12) 

CALL  GMSUBCA,BKX,AAUG,12,12) 

CALL  GMPRDCB.KU.BKU, 12,3,3) 

CALL  GMSUBCB.BKU, BAUG, 12,3) 

16  CONTINUE 

C  CALCULATE  THE  STATE  TRANSITION  MATRIX 

DO  17  1-1,12 
DO  18  J-1,12 

18  IFCI.EQ.J)  El C I , J)-DT 

17  CONTINUE 

DO  19  1-1,12 
DO  20  J-1,12 


o  o 
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20  ADT(I,J)-AAUGCI,J)*DT 
19  CONTINUE 

300  RNUM2-RNUM2+1 .0 

CALL  GMPRD(E1 ,ADT ,E2,12,12,12) 

DO  21  1-1,12 
DO  22  J-1,12 

22  E2(I,J)-E2(I,J)/'RNUM2 

21  CONTINUE 

DO  23  1-1,12 
DO  24  J-1,12 

24  IF(E2(I,J)-E1(I,J).GT. 000001)  GO  TO  25 

23  CONTINUE 
GO  TO  30 

25  DO  26  1-1,12 

DO  27  J-1,12 

27  ETAU( I , J)-ETAU( I , J)+E1 < I , J) 

26  CONTINUE 

DO  20  1-1,12 
DO  29  J-1,12 

29  E1(I,J)-E2(I,J> 

28  CONTINUE 
GO  TO  300 

30  CONTINUE 

CALL  GMPRD(AAUG,ETAU,EAT,12,12,12) 

DO  31  1-1,12 
DO  32  J-1,12 

32  IFCI.EQ.J)  EAT(I ,J)-EAT<I ,  J)+l .0 

31  CONTINUE 

CALL  GMPRD( ETAU , 8AUG , ET AUB ,12,12,3) 

PRINTOUT  OF  A,B,C,D  MATRICES 
IF  11-0  MATRICES  DILL  NOT  BE  PRINTED 
IF(Ii.EQ.0)  GO  TO  309 
URITE<5, 104) 

DO  33  1-1,12 

URITE(S, 102)  A( I , 1 ) ,A( I ,2) ,A( I ,3) ,A( I ,4) , A( I ,5) ,A( I ,6) , A( I ,7) , 
C  A( I ,B) ,A< I ,9) ,A( I , 10) , A( I , 1 1 ) ,A( I , 12) 

33  CONTINUE 

C  PRINTS  B  MATRIX 

URITE(5,105) 

DO  34  1-1,12 

34  URITE(5,103)  B(I ,1 ) ,BC I ,2) ,B( I ,3) 

C  PRINTS  C  MATRIX 

URITE(5,106) 

DO  35  1-1,10 

URITE(5,102)  C(I ,1) ,C( I ,2) ,C(l ,3) ,C(I ,4) ,C(I ,5) ,C(I ,6) 

C  ,C(I ,7) ,C(1 ,8) ,C<I ,9) ,C(I ,10) ,C(I ,1 1 ) ,C(I ,12) 

35  CONTINUE 

C  PRINTS  D  MATRIX 

URITE(5,107) 

DO  36  1-1,10 

36  UR1TE(5,103)  D( 1 ,1 ) ,D(1 ,2) ,D( 1 ,3) 

C  PRINTOUT  KX,KU  MATRICES 

URITE(5,108) 


«J  u  u  u  u  u 
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DO  37  1-1,12 

37  URITE(5,103)  KX<1,I),KX(2,I),KXC3,I) 

URITEC3.109) 

DO  38  1-1,3 

38  URITE(3,103)  KU<I ,1 ) ,KU(I ,2) ,KUCI ,3) 

C  PRINTOUT  AAUG  AND  BAUG  MATRICES 

URITE(5,118) 

DO  39  1-1,12 

URITE(5, 102)  AAUG(I ,1) ,AAUG(I ,2) ,AAUG( I ,3) ,AAUG(I ,4) , 

C  AAUG( 1,5), AAUG( 1,6), AAUG( 1,7), AAUG( 1,8), AAUG( 1,9), 

C  AAUG(I,10),AAUGCI,11),AAUGU,12> 

39  CONTINUE 
URITE<5, 111) 

DO  40  1-1,12 

40  URITE(5, 103)  BAUG( I , 1 ) ,BAUG( I ,2) ,BAUG( I ,3) 

PRINTOUT  OF  TRANSITION  MATRICIES 
URITEC5, 112)  DT, IAUG.RNUM2 
URITE(5,113) 

DO  41  1-1,12 

41  URITE(5,114)  <ETAU( I ,L) ,L>1 ,12) 

URITE(5,115> 

DO  42  1-1,12 

42  URITEC3.116)  (ETAUBU ,L),L«1 ,3) 

URITE(3,117) 

DO  43  1-1,12 

43  URITE(5, 1 14)  CEATU  ,L)  ,L-1 ,12) 

309  CONTINUE 

SET  INITIAL  VALUES 
RDCF-57 .29378 
RAD755— 75.5/RDCF 
CPYFLG-. FALSE. 

TAUR-2.0 
N2— 1 
N-0 

VTDT-0.0 
UTDT-0.0 
X(l)-.412 
X(2)-78.0 
X<  3)— 300.0 
330  CONTINUE 
T-0.0 

DO  300  1-1,4 
TGT(I)-0.0 
DTT(I)-0.0 
500  ERRT(I)-0.0 

START  THE  CLOCK  MID  READ  STICK  AND  RUDDER  VOLTAGES 
AND  CONVERT  TO  DEFLECTION  VALUES 

CALL  TSET(T) 

301  CALL  DASIN(NCHNLS,RVALS) 

CALL  TSECSCRTI) 

C  CHECK  FOR  IMPUT  FROM  KEYBOARD 

IF< .NOT. ( A8RDC( STACK, LTR)))  GO  TO  362 
IF(LTR.GE.CASE)  LTR-LTR-SHIFT 
IF(LTR.EO.'C')  GO  TO  304 
IF(LTR.EQ.'S')  GO  TO  1000 


302  CONTINUE 

U< 1 )-RVALS( 1 )*EGN 
U(2)— RVALSC2)*AGN 
UC3)-RVALS(3)*RGH 
C  CALCULATE  X  VECTOR 

DO  44  1-1,12 
EATX(I)"0.0 

44  ETAU8U(I>-0.0 
DO  45  1-1,12 

DO  46  J-1,12 

46  EATX( I )-EAT< I , J)*X( J)+EATX( I ) 

ETAUBUC I ) -ETAUBUC I )+ETAUB( 1 , 1 )*UC 1 HETAUBC I , 2)*U( 2)+ 
C  ETAUBC I ,3)*UC3) 

45  X< I )-ETAUBU< I )+EATX( I ) 

PHIREF-XC7) 

C  CALCULATE  THE  Y  VECTOR 

DO  4?  1-1,10 
CX(I)-0.0 

47  DU( I)-0.0 
DO  48  1-1,10 

DO  49  J-1,12 

49  CXCI)-CXCI)+CCI,J)*XCJ> 

48  CONTINUE 

DO  50  1-1,10 

DU< I) -D(1 ,1 >*UC 1 )+D< 1 ,2)*UC2)+DC 1 ,3)*UC3)+DUC  I ) 

50  Y( 1 )-DU< I )+CX( I ) 

C  CALCULATE  TARGET  MOTION 

DAL- 0.0 
DO  92  1-1,13 

92  DAL-DAL+AAC I )*SIN< OMEGA C I >*T+PHI T2( I ) ) 

DT2-DT/2. 

DTT( 2) — TGT<  2J/TAUR+DAL 
DTTC3)  —  .161*TGT(4) 

DTT(4)-TGT(2)+.161*TGT(3) 

J-l 

405  CONTINUE 

DO  91  1-2,4 
R1 1 -DT*DTT ( I ) 

R22-AK(J)*(R11-BKU)*ERRT(D) 

TGT(I)-TGTU)+R22 

R22-R22*3.0 

91  ERRT< I )-ERRT( I )+R22-CK( J)*R1 1 
IF( J.EQ.4)  GO  TO  406 
J-J+l 

IF( J.NE.3)  T-T+DT2 
DTT(2)— TGT(2)/TAUR+DAL 
DTT(3)  — .161*TGT(4) 

DTT<4)-TGT(2)+.161*TGT(3) 

GO  TO  405 

406  CONTINUE 
PHIREL-TGT(4) 

C  CALCULATES  TRANSFORMATION  MATRIX 

Cl-COS(X(l)) 

C2-C0S(X(6)> 

C3-C0S(X(7)) 
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c 


52 


C 


54 


C 


C4-C0S<TGT(3>) 

C5-SIN(X<1)) 

C6-SINCX<6>) 

C7-SIN<X<7>) 

C8«SIN<TGTC3>) 

C9- . 9682+ . 25*TGT ( 4) 

C10-.9682-.25*TGT<4) 

Cll-.25-.9682*TGTC4> 

C12-C2*C1*C8 

C13-C1*C6*C7*C8 

C14-C5*C6*C7 

C15-C2*C4*C7 

C16-C1*C3*C6 

C17-C2*C3*C4 

T11-C1*C2*C4+C6*C8 

T1 2* Cl 2*C9*C2*C5*C1 1 -C4*C6*C9 

T1 3-C1 2* Cl 1*C2*C5*C1 0-C4*CG*Cl 1 

T21 -C4*C7*C6*C1~C3*C5-C2*C7*C8 

T22-C13-C3*C5*C9+C14+C1*C3*C1 1+C9*C15 

T23-C13-C3*C5*C1 1-C14+C1*C3*C9+C1 1*C1 5 

T31 -C4*C1 6+C5*C7-C3*C6*C8 

T32-C8*C1S+C5*C7*C9+C14+C1*C3*C11+C9*C17 

T33-C8*C16+C5*C7*C1 1-C3*C5*C6-C1*C7*C9+C1 1*C17 

VTDT-8 . 05*T22*TGT<4)-31 . 1 8*T23*TGTC4) 

UTDT«8.05*T32*TGT<4)-31 .18*T33*TGT(4) 

X<  2) “X< 2)+VTDT 
X<  3)  -X(  3)-HJTDT 

SIGHT  POSITION  UPDATED 
C24- 300 . 8*Y ( 5 ) 

C25-300.0*Y<7) 

DO  52  1-1,8 

POSX< I ) -SX( I )+REFSX+C25 
POSYt I )-SY( I I+REFSY-C24 
CONTINUE 

CALL  GVECTU  ,POSX,POSY,DUMY,8,0000) 

TARGET  POSITION  UPDATED 
CPT-COS(PHIREL) 

SPT-SIN(PHIREL) 

C22-300.0*X(4) 

C23-300.0*X(5) 

TGTX-REFTX+C23 
TGTY- REFTY-C22 
DO  54  1-1,6 

P0SX3< I ) -TX< I )*CPT+TY< I )*SPT+TGTX 
P0SY3( I ) — TX( I )*SPT+TY< I )*CPT+TGTY 
CONTINUE 

CALL  GVECT<2,POSX3,POSY3,DUMY,6,4000) 

HORIZON  LINES  POSITIONS  UPDATED 
PHI -RAD755-PH I REF 
C20-300.0*SIN(PHI) 

C21 >300.  OnCOS( PHI ) 

DO  53  1-1,2 

POSX2( I )-»« I >*COS<PHI )+HY( I )*SINCPHI 1+510.0 
POSY2C I ) — HX< I )*SI N< PHI )+HYC I )*C0SC  PHI ) 

HX1 ( I )-P0SX2( I )+C20 
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W<2< I )-P0SX2( I ) 

HK3tI)-POSX2(l)-C20 

HY1 < I ) -P0SY2t I )+REF( 2)+C21 

HY2( I )-P0SY2( I )+REFC2) 

HY3t I )-P0SY2t I )+REF(2)-C21 
53  CONTINUE 

CALL  GVECT(3, HX1 ,HY1 ,WJMY,2 ,0000) 

COLL  GVECT  ( 4 ,  W<2 ,  HY2 ,  DUMY ,  2 , 0000 ) 

CALL  GVECT  ( 5 ,  HD<3 ,  HY3 ,  DUMY ,  2 , 0000 ) 

C  COPY  PARAMETERS  TO  DISK (EVERY  OTHER  VALUE) 

IF(ICOPY.NE.l)  GO  TO  307 
Nl— N2 

IF<N1 .GT.0)  GO  TO  308 

INDEX-INDEX+1 

XBUFt 1,1 NDEX) -T 

XBUF( 2 , I NDEX) -X( 1 )*RDCF 

XBUFC3, INDEX)-Xt2) 

XBUF(4,INDEX)-X(3) 

XBUF (5,1 NDEX ) -X 1 4 ) *RDCF 
XBUFC 6 , I NDEX) -Xt 5>*RDCF 
XBUFC7, INDEX)*X(6)*RDCF 
XBUF ( 8 , I NDEX) -X( 7)*RDCF 
XBUFO,  I  NDEX)  -X(  8)*RDCF 
XBUFt 10,1 NDEX) -X( 9)*RDCF 
XBUFt 11,1 NDEX) -Xt 1 0)*RDCF 
XBUFt 12,1 NDEX) -Xt 1 1 )*RDCF 
XBUFt 13,1 NDEX) -Xt 1 2)*RDCF 
XBUFt 14,1 NDEX) -Yt 1 )*RDCF 
XBUFt 15,1 NDEX) -Yt  2)*RDCF 
XBUF t l 6 , 1 NDEX) -Y t 3 ) *RDCF 
XBUFt 17, INDEX)-Yt4)*RDCF 
XBUFt 18,1 NDEX) -Yt  5)*RDCF 
XBUFt 1 9 , I NDEX) -Yt 6)*RDCF 
XBUFt  20, 1 NDEX) -Yt  7)*RDCF 
XBUFt  21,1 NDEX) -Yt  8)*RDCF 
XBUFt 22 , I NDEX) -Yt 9) 

XBUFt 23, I NDEX) -Ut 1 ) 

XBUFt  24 , 1 NDEX) -Ut  2) 

XBUFt25,lNDEX)-Ut3) 

XBUFt 26, I NDEX) -DAL 
XBUFt 27, I NDEX) -TGTt4)*RDCF 
XBUFt  28 , 1 NDEX) -TGT 1 2) 

IFf INDEX.LT. 7)  GO  TO  308 

CALL  FSDI Ot I WRITE , IMMED, I SECT , XBUF , I RC) 

IFt IRC.NE.0)  GO  TO  1000 
ISECT-ISECT+1 
INDEX-0 
308  CONTINUE 
N2-N1 

307  CONTINUE 

CALL  TSECStRT?) 

C  TELL  COMPUTER  TO  WAIT  TILL  REAL  TIME 

TC-RT7-RT1 

IFtTC.GT.DT)  URITEt5,119)  T,TC,DT 
CALL  TWAITtT) 
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GO  TO  301 

C********  END  OF  ITERATION  LOOP  ******** 
C******************************************** 

C  KEYBOARD  INPUT  DIECTIVES 

C  IF  KEYBOARD  INPUTS  'C',  DATA  IS  WRITTEN  TO  DISK 

304  I  COPY- 1 

CPYFLG- .TRUE. 

URITE(5,118) 

GO  TO  330 

C  STOP  SINULATION  AND  PRESENT  DESIRED  OUTPUT 

1000  CONTINUE 
CALL  GSTOP 

C  NARK  END  OF  DISK  FILE  BY  SETTING  T— 1.0 

T— 1.0 

IF(. NOT. CPYFLG)  GO  TO  2000 

INDEX-INDEX+1 

XBUFC 1 ,INDEX)-T 

CALL  FSD 1 0  C I UR I TE , I UA I T , I SECT , XBUF , I RC ) 

CALL  F*0PN<2) 

ISECT-ISECT+1 
DO  57  I SC-1, I SECT 
1 SEC- I SC-1 

CALL  FSDI 0( I READ, I WAIT, I SEC ,XBUF , I RC) 

IF( IRC.NE.0)  GO  TO  2000 
DO  58  1-1,7 

IF( I AUG.EQ.0)  GO  TO  320 
KXX( 1 )-0.0 
KXXC2)-0.0 
KXX(3)-0.0 
KUU< 1 )-0.0 
KUU(2)-0.0 
KUUC3)-0.0 
320  CONTINUE 

URITE<2,100)  (XBUF(L,I ) ,L-1 ,5) 

URITEC2, 100)  CXBUFCL, I ) ,L-6, 10) 

WRITE! 2, 100)  CXBUFCL, I ) ,L-1 1 ,15) 

WRITE! 2, 100)  (XBUF(L,I ) , L-16,20) 

WRITEC2, 100)  (XBUFCL, I) ,L-21 ,25) 

C  CALCULATES  AUGNENTED  STICK  INPUT, IF  I AUG  N.E.  0 

UAUG(  D-0.0 
UAUG(2)-0.0 
UAUG<3)«0.0 

I F ( I AUG.EQ.0)  GO  TO  322 
DO  59  J-1,3 
DO  60  K-5,13 

K»«  J) -KXX( J)+KX( J, K-l )*XBUFC K, I )/RDCF 
KXX( J)-KXX< J)+KX< J, 1 )*XBUF(2, I VRDCF+KX( J,2)*XBUF(3, I )+ 

C  KX( J,3)*XBUF(4,I) 

KUU(J)-KUU(J)+KU(J,1)*XBUFC23,I)+KU(J,2)*XBUF(24,I)+ 

C  KU( J,3)*XBUF( 25, I ) 

UAUG(J)-KXX(J)+KUU(J) 

CONTINUE 

URITE(2,100)  XBUFC 26, I ) ,XBUF(27, I ) ,UAUG( 1 ) ,UAUG(2) , UAUGC3) 
IF  (XBUFC 1,1)  .LT.  0.0)  GO  TO  1099 


60 


59 

322 


58  CONTINUE 
57  CONTINUE 
1899  URITEC2, 1 01 ) 

2800  CONTINUE 

C  THE  WRITE  STATEMENT  FORMATS 

100  F0RMAT(5X,5E15.7) 

101  FORMAT ( '  END  ',75X) 

102  FORMATC 12F9.4) 

103  F0RMAT(3F12.7> 

104  FORMAT ( ////45X ,  ’ THE  A  MATRIX V) 

105  FORMAT ( // 1 5X , ' THE  B  MATRIX'/) 

106  FORMAT ( //////////4 5X , ' THE  C  MATRIX'/) 

107  FORMAT <//15X, 'THE  D  MATRIX'/) 

108  FORMAT (///////// 10X, 'THE  KX  AUGMENTATION  MATRIX'/) 

109  FORMAT(//l 0X , 'THE  KU  AUGMENTATION  MATRIX'/) 

110  FORMAT ( //////////45X , ' THE  AAUG  MATRIX'/) 

111  FORMAT < // 1 5X , ' THE  BAUG  MATRIX'/) 

C112  FORMAT < ///// , 'THE  VALUE  OF  DT-  ' ,F5.2,/t 'THE  VALUE  OF  IAUG- 
C  C  15,/, 'THE  TOTAL  NUMBER  OF  ITERATIONS  WAS-  ',F7.1,//) 

113  FORMAT ( /////45X , ' THE  EXP  TAU  MATRIX'/) 

114  FORMAT( 12F9.4) 

115  FORMAT ( //5X , 'THE  EXP  TAU  TIMES  B  MATRIX'/) 

116  F0RMAT(3( 1X.F9.4) ) 

1 1 7  FORMAT ( //////////45X , 'THE  TRANSITION  MATRIX'/) 

118  FORMAT ( *****  DATA  DILL  BE  COPIED  TO  DISK  NOU  ****') 

119  FORMAT ( ' AT  TIME-  *,Fe.4,'  SECONDS,  THE  ITERATION  CYCLE  TIME 

C  'WAS-  ' ,F8.4,'  .WHICH  EXCEEDS  THE  SET  ITERATION  TIME 

C  'OF  ' ,F8.4, '  SECONDS',/) 

STOP 
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APPENDIX  G 
Plotting  Program 


This  appendix  contains  information  about  the  computer  routine, 
WAYPLT,  which  calculates  RMS  statistics  and  plots.  Variables  used  in 
the  program  will  first  be  defined  and  then  the  subroutines  called  will 
be  briefly  explained.  Finally,  a  computer  listing  of  WAYPLT  is  given. 


G. 1  Computer  Variables  -  small  x  represents  any  variable 


SYMBOL 

AVDESQ 


AVG 

DEF1 ,  DEF2 
DESM 

DESQ 

I  DONE 

IPLT 

IRATE 

IRMS 


DEFINITION 


Sum  of  the  variable,  squared,  and  divided  by  the 

N  2 

total  number  of  points,  ( \  x  )/N 


Average  of  the  variable. 


1 

N 

(J  x)/N 


Last  value  and  present  value  of  control  input,  u 

N  p 

Sum  of  the  variable,  7  x 

T  N  2 

Sum  of  the  variable,  squared,  |  x 

Logical  variable  to  decide  if  another  variable  should 
be  selected 

Logical  variable  to  decide  if  variables  are  plotted 
Logical  variable  to  decide  if  control  inputs  are 
calculated 

Logical  variable  to  decide  if  RMS  statistics  are 
calculated 
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Total  number  of  points,  N 

N 

RMS  value  of  the  variable,  ((£  x^/N)1^2 


Standard  deviation  of  the  variable, 

N  A  1 


((j  X)/,N)2|)1/2 


SQAVDE  Sum  of  the  variable,  divided  by  the  total  number 

N  , 

of  points,  and  then  squared,  ((^  x)/N) 

Tl,  T 2  Last  value  and  present  value  of  time,  respectively 

VNAME  Storage  array  for  variable  names  used  on  plots 

X*  Storage  array  for  plotting  points  in  the  x-direction 

XBUF  Storage  array  of  data 

Y*  Storage  array  for  plotting  points  in  the  y-direction 


♦the  integer  number  assigned  to  each  parameter  can  be  found  in  Appendix  H. 


6.2.  Subroutines 

AXES ( XSPEC , XTITLE ,  YSPEC,  YTITLE)  -  draws  and  labels  plot  axes  based  on 

scaling  factors  calculated  in  S/R  SCAN 

XSPEC,  YSPEC  -  real  number  in  form  NN.D 

if  NN  0,  NN  *  number  of  characters  in  axis  title 
if  NN  0,  no  axis  drawn 

D  ■  number  of  decimal  places  for  numbers  annotated  on  axes 
XTITLE,  YTITLE  -  array  containing  characters  of  axes  titles 

DRAW{X,  Y,  N,  IDS)  -  draws  a  line  between  a  sequential  set  of  coordinates, 

or  draws  a  point  at  each  coordinate 

X,Y  -  real  or  integer  arrays  of  coordinates 
N-  integer  value,  whose  absolute  value  equals  the  number  of 
coordinates  if  N  *  0,  plotting  is  terminated 
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IDS  -  integer  data  specification,  in  the  form  txyc 
t  «  line  type 

t  *  0,  draws  straight  lines  between  points 
t  *  9,  draws  straight  line  to  the  last  point 
x,y  -  type  of  data  in  X,Y  arrays 

0  =  real  value  that  does  not  have  to  be  scaled 
4  =  real  value  that  must  be  scaled 
C-line  width 

0  =  previous  width 

F$0PN{1)  -  opens  up  the  assigned  storage  location  so  the  information 
can  be  manipulated  by  FORTRAN  programming,  such  as  READ  and  WRITE  commands 
FREMAT  -  allows  usage  of  FORTRAN  free-format  READ  and  WRITE  statements 
LINK  -  terminates  plotting  program  and  allows  computer  system  to  call 
VPLOT  program,  which  generates  the  plots 

SCAN(X,  Y,  N,  IDS)  -  scans  the  data  in  the  storage  arrays  and  calculates 

scaling  factors  to  map  coordinates  onto  the  electrostatic  printer/plotter 

X,  Y  -  real  or  integer  arrays  of  coordinates 
N  -  integer  value,  whose  absolute  value  equals  the  number  of 
coordinates  if  N  <0,  determines  minimum  and  maximum  values 
of  coordinates  and  calculates  scale  factors 
IDS- integer  data  specification,  in  the  form  xyO 
x,y  -  type  of  data  in  X,Y  arrays 
4  =  real  value  that  must  be  scaled 


100 


C**  UAYPLT  -  PLOTS  DISPLY  DATA  VS  TIME  F 

DIMENSION  X<900) ,Y(900) ,XBUF(30) ,VNAME<33) 

DOUBLE  PRECISION  DESM , DESQ , flVDESO , SQAVDE , SDDE , RMDE , AVG 
DATA  VNAME/'TIME' ,'PSI  VT  ','  UT  ','T  EL', 

C  'T  AZ'.'TETA  ','PHI  ','ALPH','  Q 

C  'BETA','  P  ','  R  ','E  EL'.'E  RT' , 

C  'E  AZ'.'EA  R'.'S  EL','SE  R','S  AZ', 

C  'SA  R'.'TACC','  DE  DA  ','  DR  ', 

C  'DAL  '.'PHIT'.'DE  A' ,'DA  A'.'DR  A', 

C  'DE  R'.'DA  R' ,'DR  R'/ 

5  URITE(7,10) 

10  FORMAT ( 

C 
C 
C 
C 
C 
C 

CWWtWWt  VALUES  OF  IY:  1-T,  2-PSI ,  3-VT,  4-UT,  5-TGT  EL  ******* 
cmtstnst  6s AZ  T,  7-  THETA,  B-PHI,  9-ALPHA,  10-Q,  11 -BETA  ******* 


C*»»*****»*  12-P,  13-R,  14-EL  ERR,  15-EL  ERR  RT,  16-E  AZ  ******* 
C***i**»«»«*  17- AZ  ERR  RT,  18-SGT  EL  ERR,  19-SGT  E  RT  ******* 
CSSS*******  20-S  AZ  ERR,  21-SAZ  E  RT,  22-TGT  ACC,  23- DE  ******* 
C**********  24-DA,  25-DR,  26-DA  TGT  ,27-PHI  TGT  ******* 
C**»**»»***  28- DE  AUG,  29- DA  AUG,  30- DR  AUG  ******* 
C**********  31 -DE  RATE,  32- DA  RATE,  33- DR  RATE  ******* 


P#  L  FREMAT 

READ( 7,11)  IX, IY.IPLT.IRMS, IRATE 

11  FORMAT<V) 

ICNT-0 

IY2-IY 

IF( IRATE. EO.l)  IY2-IY+8 

Tl-0.0 

DEF1-0.0 

DESM- 0.0 

DESQ- 0.0 

CALL  FSOPN(l) 

12  READ(1 ,20)  (XBUF(L) ,L-1 ,5) 

READ( 1 ,20)  (XBUF(L) ,L-6, 10) 
READ(1 ,20)  (XBUF(L) ,L-1 1 ,15) 
READd ,20)  (XBUF(L) ,L-16,20) 
READd ,20)  (XBUF(L) ,L-21 ,25) 
READd, 20)  (XBUF(L)  ,L-26,30) 
IF(XBUFd).EQ.-l)  GO  TO  100 
ICHT-ICNT+1 

IFdRATE.EO.0)  GO  TO  15 
X<ICNT)-XBUF(IX) 

T2-XBUFCIX) 

DT-T2-T1 

T1-T2 

DEF2-XBUFC IY) 

DELTA- DEF2-DEF1 
DEF1-DEF2 
YCICHT)-DELTA/DT 
60  TO  16 


Wow  ENTER  IX, IY.IPLT.IRMS, IRATE  wow',/ 

Wow  IX, IY-SUBSCRIPT  OF  PARA  BEING  PLOTTED  wow',/ 
wow  IPLT-0  IF  NO  PLOTS,  IPLT-1  IF  PLOT  REO  wow',/ 
wow  I RMS-0  IF  RMS  NOT  CALC,  IRMS-1  IF  CALC  Wow' ,/ 
wow  I RATE- 1  IF  STICK  RATES  CALC, ELSE  I RATE-0  wow',/ 

WuWWnWWoWWoWWtWOWWnWWuWWcWtWoWWoWWoWWoW*') 


11/24/B2 
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15  CONTINUE 

X( ICNT)«X8UF( IX) 

Y<ICNT)-X8UF(IY) 

16  CONTINUE 

IF( IRMS.EQ.0)  GO  TO  30 
DESM-DESmY(ICNT) 

DESQ- DESO+Y ( ICNT)*Y( ICNT) 

30  CONTINUE 

20  F0RMAT(5X,5E15.7) 

GO  TO  12 

C  CALCULATE  RMS  AND  STD  DEV,  IF  REQ 

100  IFURMS.EQ.0)  GO  TO  40 
PTS-FLOATCICNT) 

AVG-DESrVPTS 
AVDESQ-DESQ/PTS 
SQAVDE-CDESM/PTS)**2.0 
SDDE-SQRTC ABS( AVDESO-SQAVDE ) ) 

RMDE-SQRT ( AVDESQ) 

URITEC5.21 )  VNAMEC IY2) ,AVG,VNAME( IY2) , RMDE, VNAMEC I Y2) ,SDDE,ICNT 

21  FORMAT  (  '  THE  AVG  VALUE  OF  *,A4,'  IS-' ,F1 1 .5,/', 

C  '  THE  RMS  VALUE  OF  ',A4,'  IS-' ,F1 1 .5,/, 

C  '  THE  STD  DEV  OF  *,A4,'  IS-',F11.5,/ 

C  '  THE  NUM  OF  POINTS-' , 15,/) 

40  CONTINUE 

C  START  PLOTTING 

IFUPLT.EQ.0)  GO  TO  41 
CALL  SCAN<X,Y, -ICNT, 0440) 

CALL  AXES<4 .4 ,VNAME( IX) ,4 .4 , VNAMEC I Y2) ) 

CALL  DRAUCX,Y, ICNT, 0440) 

CALL  DRALKO. 0,0. 0,1, 9000) 

41  CONTINUE 
URITEC7.13) 

13  FORMAT C '************************«**iMi*****»iM<*«***********  , 

C  *****  ENTER  IDONE  ****',/, 

C  *****  IF  IDONE-0, START  PLOTTING  ****’,/, 

C  *****  IF  I DONE- 1, SELECT  NEU  PARAMETER  TO  PLOT  ****',/, 

C  '**************************************************'  t/) 

RE ADC  7, 14)  I DONE 

14  FORMATCV) 

IFC I  DONE. EQ. 1 )  GO  TO  5 
CALL  LINKC'VPLOT  *) 

STOP 

END 
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APPENDIX  H 
Simulation  Results 


Appendix  H  contains  the  identification  procedure  used  to  store  the 
simulation  data  and  defines  the  variables  stored.  In  addition,  the 
complete  set  or  RMS  statistics  used  to  analyze  the  simulation  are 
tabulated.  Finally,  the  comments  of  each  individual  Mho  flew  the 
simulation  are  given. 

HI.  Datafile  Identification  Procedure 

The  four  individuals  who  took  part  in  the  simulation  are  identified 
by: 

(1)  DJB  -  Lt.  Colonel  in  the  US  Air  Force,  with  4000  hours  of 
flying  time  in  helicopters,  transports,  and  as  a  test  pilot. 

(2)  DAW  -  Captain  in  the  US  Air  Force,  with  1200  hours  of  KC-135 
flying  time. 

(3)  BJH  -  Purdue  graduate  student  with  25  hours  of  civilian 
flying. 

(4)  WAY  -  Captain  In  the  US  Air  Force,  with  900  hours  of  F-4 
flying  time. 

The  data  obtained  In  the  simulation  was  stored  on  disk  using  the 
storage  location  identifier  of  XXXA8C,  where 

XXX  *  Individuals  Initials,  given  above 
A  *  augmentation  Identifier 
If  A  «  0,  no  augmentation 
if  A  ■  1,  augnentatlon 


-  4.  f. 


B  *  set  equal  to  zero  for  this  simulation 
C  *  the  simulation  run  number  for  that  particular  augmentation 
configuration 

For  example,  DJB101  contains  data  for  Individual  DJB's  first  simulation 
run,  with  augmentation.  The  variables  stored  in  the  datafiles  are  given 


below 


Storage 

Location 

Analytic 

Identifier 

Computer 
Variable  Name 

Meaning 

Units 

1 

t 

T 

Time 

sec 

2 

Aip 

X(l) 

Heading  difference 

deg 

3 

VT 

X(2) 

Target  velocity  component 

ft/sec 

4 

WT 

X(3) 

Target  velocity  component 

ft/sec 

5 

eEL 

X(4) 

Target  elevation  angle 

deg 

6 

*AZ 

X(  5) 

Target  azimuth  angle 

deg 

7 

e 

X(6) 

Attacker  pitch  angle 

deg 

8 

4> 

X(7) 

Attacker  bank  angle 

deg 

9 

a 

X(8) 

Attacker  angle  of  attack 

deg 

10 

q 

X(9) 

Stability  axis  angular  rate 

deg/sec 

11 

B 

X(10) 

Attacker  sideslip  angle 

deg 

12 

p 

X(ll) 

Stability  axis  angular  rate 

deg/sec 

13 

r 

X(  12 ) 

Stability  axis  angular  rate 

deg/sec 

14 

eEL 

Y(l) 

Elevation  error 

deg 

15 

eEL 

Y(2) 

Elevation  error  rate 

deg/sec 

16 

eAZ 

Y(3) 

Azimuth  error 

deg 

17 

*AZ 

Y(4) 

Azimuth  error  rate 

deg/sec 

18 

XEL 

Y(5) 

Sight  elevation  angle 

deg 

19 

*EL 

Y(6) 

Sight  elevation  angle  rate 

deg/sec 

20 

XAZ 

Y(7) 

Sight  azimuth  angle 

deg 

Storage 

Location 

Analytic 

Identifier 

Computer 
Variable  Name 

Meaning 

Units 

21 

*AZ 

Y(8) 

Sight  azimuth  error  rate 

deg/sec 

22 

Ay 

Y(9) 

Attacker's  lateral  accel. 

ft/sec 

23 

6e 

at 

u(D 

Elevator  deflection 

In. 

24 

*st 

u(2) 

Aileron  deflection 

In. 

25 

6Vd 

u(3) 

Rudder  deflection 

In. 

26 

\ 

DAL 

Target  aileron  deflection 

radl an s 

27 

TGT(4) 

Target  perturbation  bank 
angle 

deg 

28 

W'* 

Augmentation  elevator 
deflection 

In. 

29 

Maug 

U.ug<2) 

Augmentation  aileron 
deflection 

In. 

30 

H2.  RMS_ 

\ug  i 

1 

Statistics 

W3) 

Augmentation  rudder 
deflection 

In. 

This  section  contains  tabulated  results  of  RMS  statistics  calculated 


for  each  Individual. 
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Table  18 

DJB  Simulation  Results  (RMS) 

PARAMETER  DJB001  DJB002  DJB101  DJB 102 


eEL 

.79 

.59 

.31 

.33 

CAZ 

3.55 

2.98 

1.37 

1.31 

VT 

19.52 

16.51 

6.17 

6.11 

«T 

28.64 

15.71 

13.77 

27.76 

a 

.79 

.49 

.29 

.67 

B 

.32 

.28 

.35 

.40 

XEL 

.84 

.42 

.30 

.84 

XAZ 

.54 

.54 

•17 

.20 

.1 

0 

8.71 

10.96 

10.96 

9.11 

$ 

10.56 

11.48 

10.66 

9.54 

P 

11  66 

10.29 

3.04 

2.66 

q 

1.78 

1.28 

.57 

1.10 

r 

1.04 

1.18 

.47 

.44 

6E 

.22 

.13 

.10 

.23 

6a 

.68 

.58 

.18 

.18 

6r 

.36 

.19 

.36 

.36 

!e 

.79 

.57 

.36 

.41 

*A 

2.60 

2.58 

.55 

.48 

1.87 

1.01 

.96 

.92 

Ay 

3.95 

2.69 

2.36 

1.61 

♦T( AUG/RHS) 

-.39/5.28 

-.40/5.23 

-.40/5.23 

-.40/5.21 

i 


’*•*»**£<( 


21.63 

32.23 


17.40 

19.35 


10.63 


10.94 

11.60 


*t(AUG/RMS)  -.40/5.22  -.40/5.26 


--■ire**  HHECi*  ~  -  f&*r‘ i 


Table  21 

WAY  Simulation  Results  (RMS) 


i 


PARAMETER 

WAY001 

WAY002 

WAY101 

WAY102 

eEL 

1.01 

.98 

.39 

.51 

CAZ 

4.41 

4.05 

1.24 

1.11 

VT 

21.52 

19.64 

6.17 

5.57 

“t 

23.59 

26.05 

23.24 

28.56 

a 

.80 

!  -81 

.49 

.69 

6 

•  71 

t 

.67 

.45 

.40 

XEL 

.60 

.71 

.69 

.86 

XAZ  | 

.60 

.57 

.16 

.15 

0  i 

10.36 

9.55 

10.10 

9.10 

4> 

i 

| 

:  11.85 

10.80 

10.00 

9.60 

i 

P  i 

13.82 

13.57 

2.65 

2.03 

q 

2.79 

2.52 

.88 

1.12 

r 

2.86 

2.80 

.30 

.24 

6E 

00 

CM 

• 

.27 

.17 

.25 

6a 

.74 

.70 

.13 

.21 

“0 

‘0 

.42 

.37 

1.64 

1.55 

.20 

.29 

5* 

3.41 

3.44 

.29 

.31 

-o 

'0 

.47 

.26 

Ay 

3.37 

3.15 

3.09 

3.41 

♦t(AU6/RMS) 

-.40/5.22 

-.31/5.09 

-.40/5.22 

-.40/5.21 

40/5.26 
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H.3  Pilot  Comments 

This  section  contains  the  comments  of  each  subject. 

Pilot  DJB: 

Was  the  augmentation  system  helpful?  Yes. 

Any  undesireable  or  annoying  characteristics?  No 
How  would  you  rate  the  task  difficulty,  from  1-10,  where 
1  *  no  inputs  required 
5  =  reasonable  difficulty  for  the  task 
10  =  uncontrollable 
Unaugmented  system?  5-6 
Augmented  system?  2-3 
Comments:  Fairly  realistic  task 

Pilot  8JH: 

Was  the  augmentation  system  helpful?  Yes 

Any  undesireable  or  annoying  characteristics?  No 
How  would  you  rate  the  task  difficulty,  from  1-10,  where 
1  =  no  inputs  required 
5  =  reasonable  difficulty  for  the  task 
10  *  uncontrollable 
Unaugmented  system?  6-7 
Augmented  system?  2-3 

Comments:  Fun  simulation.  Without  augmentation,  you  had  to 
anticipate  the  target  motion  by  watching  its  bank. 


no 


Pilot  DAM: 

Was  the  augmentation  system  helpful?  Yes^ 

Any  undesireable  or  annoying  characteristics?  No 
How  would  you  rate  the  task  difficulty,  from  1-10  where 
1  =  no  inputs  required 
5  =  reasonable  difficulty  for  task 
10  =  uncontrollable 
Unaugmented  system?  7 
Augmented  system?  3 

Comments:  Initially,  when  using  the  rudder  and  aileron  to  track  the 
target,  the  two  control  surfaces  were  actually  fighting 
each  other. 


Pilot  WAY: 

Was  the  augmentation  system  helpful?  Yes_ 

Any  undesireable  or  annoying  characteristics?  No 
How  would  you  rate  the  task  difficulty,  from  1-10,  where 
1  =  no  inputs  required 
5  =  reasonable  difficulty  for  the  task 
10  =  uncontrollable 
Unaugmented  system?  7 
Augmented  system?  2 

Comments:  Difficult  to  maintain  azimuth  control  without  augmentation. 


